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The recent observation of vr vr at the B factories with a surprisingly large 



branching ratio represents a challenge for theory, and complements the amazingly 
small Bd — > vr^vr" rate. We point out that all puzzling B vrvr features can be ac- 
• commodated in the Standard Model (SM) through non-factorizable hadronic inter- 

ference effects, extract the relevant parameters, and predict the CP asymmetries of 
Bd TT^n^. Using then SU{3) flavour-symmetry and plausible dynamical assump- 
tions, we fix the hadronic B — > ttK parameters through their B vrvr counterparts, 
5^ I and determine the CKM angle 7, with a result in remarkable accordance with the 

usual fits for the unitarity triangle. We may then analyse the B irK system in 
the SM, where we find agreement with the experimental picture, with the exception 
of those observables that are significantly affected by electroweak (EW) penguins, 
thereby suggesting new physics (NP) in this sector. Indeed, a moderate enhance- 
ment of these topologies and a large CP-violating NP phase allow us to describe any 
currently observed feature of the B — > ttK modes, and to predict the CP-violating 
Bd tt^Kq, observables. If we then restrict ourselves to a specific scenario where 
NP enters only through penguins, we obtain a link to rare K and B decays, 
where the most spectacular NP effects are an enhancement of the Kj^ n^vi? 
rate by one order of magnitude with BR(i^rL tt^vv) ~ 4BR(i^+ — > vr+z^P), 
BR(i^L TT^e+e") = 0(10^^°), (sin 2/3)7r;/p < 0, and a large forward-backward 
CP asymmetry in Bd K*fj,~^fi~. We address also e'/e and other prominent 
decays, including B (f)K and B J/ipK modes. 
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1 Introduction 



In this decade, dedicated B- and i^-decay experiments aim at stringent tests of the 
flavour dynamics of the Standard Model (SM) and in particular of the Kobayashi- 
Maskawa mechanism of CP violation [1]. The central target of these studies, which 
will hopefully shed light on new physics (NP), is the well-known unitarity triangle (UT) 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, with its three angles a, P and 7 
(for detailed reviews, see [2, 3]). Thanks to the SLAC and KEK B factories with their 
detectors BaBar and Belle, respectively, mixing-induced CP violation is now an estab- 
lished effect in the i?-meson system. The corresponding determination of sin (p^ through 
the "golden" mode B^ —>■ J/ipK^, where (pd denotes the B^-B^ mixing phase {(pd = 2j3 
in the SM) , agrees remarkably well with the CKM fits [4] that include in particular the 
size of the well-known "indirect" CP violation in Ki^ — > tttt decays. In spite of this 
tremendous success of the SM, one should realize that only a handful of CP-violating 
decays and rare B and K decays have been measured and it is to be seen whether some 
modifications of the SM picture of flavour dynamics and CP violation will be required 
in the future when the data improve. For instance, the BaBar and Belle data on 0^ 
still admit two solutions with 0^ ~ 47°, in accordance with the SM, and 133°, pointing 
towards NP contributions. This ambiguity can be resolved by measuring the the sign of 
cos0d [5], which is in progress at the B factories, even though it is a challenging task [6]. 

While testing the SM and its possible extensions it is essential to consider simultane- 
ously as many processes as possible. Only in this manner can the parameters of a given 
theory be fully determined; having them at hand, predictions for other observables can 
be made. In this enterprise correlations between various observables play an important 
role, as they may exclude or pinpoint a given extension of the SM even without a detailed 
knowledge of the parameters specific to this theory. 

Interestingly, the current S-factory data for a number of processes indicate potential 
inconsistencies with the SM description of CP violation and flavour dynamics that may 
suggest the presence of NP contributions and/or the deficiencies in our understanding 
of hadron dynamics that necessarily enters the analyses of non-leptonic B decays. In 
particular: 

• The BaBar and Belle collaborations have very recently reported the observation of 
Bd — > tt^tt" decays with CP-averaged branching ratios of (2.1±0.6±0.3) x 10^^ and 
(1.7 ± 0.6 ± 0.2) X 10~^, respectively [7, 8]. These measurements represent quite 
a challenge for theory. For example, in a recent state-of-the-art calculation [9] 
within QCD factorization [10], a branching ratio that is about six times smaller is 
favoured, whereas the calculation of B^ — > 7r~^7i~ points towards a branching ratio 
about two times larger than the current experimental average. On the other hand, 
the calculation of S+ — > 7r+7r° reproduces the data rather well. This "B — > tttt 
puzzle" is reflected by the following quantities: 




2 I BR{B+ tt+tt") + BR{B- 

[bR(SO ^ TT+TT-) + BR(SO 



TT+TT )J Tb+ 



TT 7r°)1 T-fiO 



= 2.12 ±0.37 



(1.1) 
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' BR(E° ^ ttV) + BR(^° ^ 7r°7r°) ' 

_BR(SO ^ TT+TT-) + BR(SO ^ TT+TT-)_ 



- 0.83 ±0.23. 



(1.2) 
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In order to calculate the numerical values, we have used tb+ /t^o = 1.086 ± 0.017 
and the most recent compilation of the Heavy Flavour Averaging Group (HFAG) 
[11], adding the errors in quadrature. The central values calculated within QCD 

factorization [9] give R'^_ = 1.24 and K^q = 0.07. As was recently pointed out 
[12], these data indicate important non-factorizable contributions rather than NP 
contributions, and can be perfectly accommodated in the SM. 



In the B ttK system, the CLEO, BaBar and Belle collaborations have measured 
the following ratios of CP-averaged branching ratios [13]: 



Rr 



BR{D- 



Vv'-) + BR(Z?- 



7r+X0) + BR(S- ^TT-XO). 



BR{B+ 



1.17 ±0.12 (1.3) 



= 0.76 ±0.10, (1.4) 



with numerical values following from [11]. As noted in [14], the pattern of i?c > 1 
and - in particular - R^ < 1, which is now consistently favoured by the separate 
BaBar, Belle and CLEO data, is actually very puzzhng. On the other hand, the 
quantity [15] 



R 



BR(E° ^ 7r-K+) + BR(E° ^ tt+R-) 
BR{B+ TT+i^o) + BR{B- 7r-i?o) 



rB+_ 



0.91 ±0.07 



:i.5) 



does not show any anomalous behaviour. Since Rc and 7?„ arc affected significantly 
by colour-allowed electrowcak (EW) penguins, whereas these topologies may only 
contribute to R in colour-suppressed form, this "B — > ttK puzzle" may be a mani- 
festation of NP in the EW penguin sector [12, 14, 16], offering an attractive avenue 
for physics beyond the SM to enter the B — > nK system [17]-[20]. 

Another potential discrepancy with the SM expectation is indicated by the decay 
Bd (l)Ks. Within the SM, this transition is governed by QCD penguins [21] and 
receives sizeable EW penguin contributions [22, 23], so that it may well be affected 
by NP. The current experimental status of the CP- violating B^ (pKa, observables 
is given as follows [24, 25]: 



^cp( 



Ba 



cj>Ks 



-0.38 ± 0.37 ± 0.12 (BaBar) 
+0.15 ± 0.29 ± 0.07 (Belle) 

-0.45 ± 0.43 ± 0.07 (BaBar) 
+0.96 ± 0.50l^:J^ (Belle), 



;i.6) 



:i.7) 



where we have employed the same notation for the direct and mixing-induced CP 
asymmetries Acp{Bd — > (pKg) and AQp{Bd 4>Ks), respectively, as in [2].^ In 
the SM, the direct asymmetry is O(A^), i.e. tiny, whereas the mixing-induced CP 



^This notation may differ in signs from the one used by some authors. However, the exphcit example 
in (3.22) and (3.23) should allow for an easy comparison of signs. 
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asymmetry is - up to corrections of O(A^) - equal to the one of the "golden" decay 
Bd — J/i'Ko, [26]-[30], which is measured to be 

A'^^iBa J/^Ks) = -0.736 ± 0.049. (1.8) 

Consequently, we may well arrive at a discrepancy with the SM description of CP 
violation, although the experimental situation is of course very unclear at present. 

In view of significant experimental uncertainties, none of these exciting results is con- 
clusive at the moment, but it is legitimate and interesting to take them seriously and to 
search for possible origins of these "signals" for deviations from the SM expectations. As 
we are dealing here with non-leptonic decays, the natural question arises of whether these 
signals originate in the NP contributions or/and result from our insufficient understand- 
ing of the hadron dynamics. The purpose of the present paper is to develop a strategy 
that would allow us to address this question in a systematic manner once the experi- 
mental data on the relevant non-leptonic B decays and rare B and K decays improve. 
In order to illustrate this strategy in explicit terms, we shall consider a simple extension 
of the SM in which NP enters dominantly through enhanced penguins involving a 
CP-violating weak phase. As we will see below, this choice is dictated by the pattern of 
the data on the B nK observables and the great predictivity of this scenario. It was 
first considered in [31]-[33] to study correlations between rare K decays and the ratio 
e'/e measuring direct CP violation in the neutral kaon system, and was generalized to 
rare B decays in [34]. Here we extend these considerations to non-leptonic 5- meson de- 
cays, which allows us to confront this extension of the SM with many more experimental 
results. Our strategy consists of three interrelated steps, and has the following logical 
structure: 

Step 1: 

Since B — > tttt decays and the usual analysis of the UT are only insignificantly affected by 
EW penguins, the B tttt system can be described as in the SM. Employing the SU{2) 
isospin flavour symmetry of strong interactions and the information on 7 from the UT 
fits, we may extract the relevant hadronic parameters, and find large non-factorizable 
contributions, which are in particular refiected by large CP-conserving strong phases. 
Having these parameters at hand, we may then also predict the direct and mixing- 
induced CP asymmetries of the B^ tt^tt^ channel. A future measurement of one of 
these observables allows a determination of 7. 

Step 2: 

If we use the SU{3) flavour symmetry and plausible dynamical assumptions, we may 
determine the hadronic B — > nK parameters through the B — > tttt analysis, and may 
calculate the B nK observables in the SM. Interestingly, we find agreement with the 
pattern of the S-factory data for those observables where EW penguins play only a minor 
role. On the other hand, the observables receiving significant EW penguin contributions 
do not agree with the experimental picture, thereby suggesting NP in the EW penguin 
sector. Indeed, a detailed analysis shows that we may describe all the currently available 
data through sizeably enhanced EW penguins with a large CP- violating NP phase around 
—90°, in the spirit of the NP scenario considered here. A crucial future test of this 
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scenario will be provided by the CP-violating Bd — > ir^Ks observables, which we may 
predict. Moreover, we may obtain valuable insights into SU (3)-breaking effects, which 
support our working assumptions, and may also determine the UT angle 7, with a result 
in remarkable agreement with the well-known UT fits. 

Step 3: 

In turn, the sizeably enhanced EW penguins with their large CP-violating NP phase 
have important implications for rare K and B decays as well as e'/e, where the new 
weak phase plays a particularly important role. Interestingly, several predictions differ 
significantly from the SM expectations and should easily be identified once the data 
improve. Similarly, we may explore specific NP patterns in other non-leptonic B decays 
such as Ba 4>Ks. 

The most interesting results of this study have recently been summarized in [12]. Here 
we discuss the details of our analysis, present several additional results, and propose 
other methods that will be useful for the confrontation of the forthcoming data with 
the SM and the search for possible indications of NP. The outline of this paper is as 
follows: in Section 2, we discuss in detail our scenario of NP, using low-energy effective 
Hamiltonians as the starting point. In the subsequent three sections, we execute the 
three steps described above. In Section 3, we discuss the S — > tttt system in detail. We 
then move on to the B — > ttK system in Section 4, and explore the impact of this study 
on rare K and B decays and e'/e, as well as the implications of the current data on rare 
decays on the B — > 71 K observables, in Section 5. As the last element of our analysis, 
we investigate in Section 6 the effects of our NP scenario on the prominent S-meson 
decays of the kind B — > (pK and B J/ipK. Finally, we summarize our conclusions in 
Section 7. A compendium of the most relevant formulae for B tttt, Bt^^-^ ttK and 
Bs K^K^ observables and some technical details can be found in Appendices A-D. 

Our paper discusses a large number of observables and it is useful already now to 
list them, as well as the parameters they depend on, and indicate where the explicit 
expressions for them can be found in our paper. We do this in Table 1, where we 
also present our predictions (TH), list the input values and the present experimental 
values (EXP) for the observables in question. It should be emphasized that, within our 
approach, the hadronic B(^s) ~^ ^ri^ and Bg — > K^K~ parameters can be calculated in 
terms of the parameters of the S — > tttt system. The relevant formulae are given in 
(4.17), (4.18), (4.20) and (4.21). A compendium of the formulae for all observables of 
Table 1 is given in Appendix A. 

2 A Simple Scenario for New Physics 
2.1 General Structure 

The scenario of NP with the dominant Z°-penguin contributions presented here was 
first considered in [31]-[33], where correlations between rare K decays and e'/e were 
studied. It was generalized to rare B decays in [34]. The new feature of our analysis is 
a simultaneous study of rare K and B decays, in addition to non-leptonic B decays, in 
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Quantity 


Eq. 


Parameters 


TH 


EXP 




(3.20) 


X, A, d, 6 


input 


2.12 ± 0.37 


-"-00 




rr- \ rl 

X, LI, a, 17 


input 


U.oo It U.zo 




(3.25) 


d, e 


input 


-0 38 ± 16 




(3.26) 


d, 6, (pd 


input 


0.58 ±0.20 




(3.69) 


X, A, d, 9 


-0.4i;0;?^ 






(3.70) 


X A d 6 6^ 


-0 55+^1^ 




it 


(A '^0\ 


r, 5 


n QAq+o oas 

u.y'±o_Q 


n Q1 -1- n 07 

u.yi zt u.u / 




t A QR") 
1^4. ou ) 




i.UU_o.o8 


1 1 7 -1- n 1 

i.i / zt U.iZ 


p 




0, Vc, Oc, q, (JJ, (p, Pn, 


rv 09+0.12 


n 7fi -1- n in 






/ . 




n nor, 1 A 098 
u.uyj ± u.uzo 








~ 


-0 02 ± 06 


A^^{B^^7r^K^) 


(4.50) 


Tr.. Sr.. a. cu. (b 


0.03+no^ 


0.00 ± 0.07 


A'^l>{B,^7r''Ks) 


(4.59) 


r, S, rc,Sc,q,uj,(l),Pn,en 






A'SpiB.^n^Ks) 


(4.60) 


r, S,rc,Sc,q,uj,(f),pn,9n,(pd 


-o.98;h^ 


-0.48l°;l« 


A''Sl>{Bs^K+K-) 


(4.70) 


d, 9 






A^f{B,^K+K-) 


(4.71) 


d, 9, (j)s 


-o.istH? 






(4.83) 


r, S 


-0.381°:}^ 





Table 1: A guide to master formulae. We suppress the dependence on 7 as it enters all 
quantities. For those quantities that depend on the EW penguin parameters q and cp, the 
rare decays constraints as discussed in Section 5 have been employed in the calculation of our 
theoretical values (TH). 

particular B — > tttt, ttK modes. 

It should be noted that in [31]-[34] model-independent analyses and studies within 
particular supersymmetric scenarios were presented. As we will see below, in the present 
analysis, we determine the size of the enhancement of the Z'^-penguin function C and the 
magnitude of its complex phase by the B ttK data and, consequently, our predictions 
are more specific than was possible in these papers. In what follows, we will recall the 
basic ingredients of this scenario in a notation suitable to our analysis. 

In order to discuss the non-leptonic i?-meson decays and rare K and B decays in 
a particular extension of the SM, and to investigate correlations between the NP con- 
tributions to different observables, it is essential to formulate the theory with the help 
of an effective Hamiltonian. The effective weak Hamiltonian for AF — 1 decays with 
F = S, B is generally given as follows [35] : 

"^eflf = ^ X] VcKMCi{p)Qi, (2.1) 

where Gp is the Fermi constant and V^km the relevant CKM factors. Next, the Qi 
are local operators, which govern the decays in question, and the Ci{p) are the corre- 
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spending Wilson coefficients, which describe the strength with which a given operator 
enters the Hamiltonian. The latter coefficients can be calculated in renormalization- 
group-improved perturbation theory and carry, in particular, the information about the 
physics contributions at scales higher than ji, which is usually chosen to be 0{mi)) and 
(9 ((1-2) GeV) for B decays and K decays, respectively. Thus, the Ci include the top- 
quark contributions and contributions from other heavy particles if extensions of the SM 
are considered. Consequently, the Ci(/i) generally depend on the top-quark mass and 
also on the masses and couplings of the new particles. This dependence can be found 
by evaluating box and penguin diagrams with W-, Z-, top- and new-particle exchanges 
and properly including short-distance QCD effects. 

The amplitude for the decay of a given meson M = K,B, ... into a final state F = 
TTTT, nK, nuu, ... is then simply given by 

A{M ^ F) = (F|Kff|M) = ^J2^^KMCi{fi){Fmf^)\M), (2.2) 

where the {F\Qi{fj,)\M) are the hadronic matrix elements of the operators Qi between 
M and F, evaluated at the renormalization scale /i. As demonstrated in [36, 37], the 
formula (2.2) can be cast into the following useful expression: 

>l(Decay) = Po(Decay) + ^ P^(Decay) Fr{v). (2.3) 

r 

To this end, we choose n = no = 0{M\v ^mt) in (2.2) and rewrite the Ci{fio) as linear 
combinations of the so-called master functions Fr{v), which result from various pen- 
guin and box diagrams with heavy particle exchanges and v denoting collectively the 
parameters of a given model. In the SM, the functions Fr{v) reduce to the Inami-Lim 
functions [38], with v — rn^/M^. The term Pq summarizes the contributions stemming 
from light internal quarks, such as the charm and up quarks, and the sum incorporates 
the remaining contributions. The general properties of the Pq, Pr and Fr have recently 
been discussed at length in [39]. We recall here only the following features: 

• Fr{v) are process-independent, universal functions that depend on the particular 
model considered. NP enters the decay amplitudes only through these functions. 

• Pq and Pr are process- dependent quantities. In particular, they depend on the 
hadronic matrix elements of the operators Qi.^ 

In the models with "minimal flavour violation" (MFV), as defined in [39, 40], the set of 
the Fr{v) consists of seven functions 

S{v), X{v), Y{v), Z{v), E{v), D'{v), E'{v), (2.4) 

which arc discussed in detail in [39]. The important property of the functions in (2.4) 
is that they are real-valued, so that the CP-violating effects are governed entirely by 
the complex CKM phase hiding in the parameters Pr- Other definitions of MFV can be 
found in [41, 42]. 

^ These quantities should not be confused with the penguin parameters in Sections 3 and 4. 
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2.2 Going Beyond MFV 

In the present paper, we would like to investigate a scenario along the lines of [31]-[34], 
which goes beyond the usual MFV framework. In particular, we consider the presence of 
a single weak complex phase in the functions Fr{v), leading to additional CP- violating 
effects that are not present in models with MFV. In order to make the discussion as 
simple as possible, we will choose the functions S, E, D' and E' to be real-valued, as 
in the MFV case. Moreover, we will assume that the functions -E, D' and E' are not 
affected significantly by NP contributions, so that D' and E' , which result from magnetic 
penguin diagrams, describe with their SM values automatically the B — > Xs^ decay. The 
function E resulting from QCD penguins and entering the Wilson coefficients of penguin 
operators at scales 0{M\y) is phenomenologically not very relevant and can safely be 
set to its SM value. The point is that the contributions of QCD penguins to essentially 
all decays are dominated by renormalization-group effects for scales below Mw and are 
included in the coefficient Pq in (2.3). 

Concerning the AF = 2 box function S{v), it enters only the CP violation parameter 
Ek and the B^^-B^^ mass differences AMd^s- Being real, it does not introduce any new 
complex phases in the standard analysis of the UT that is based on \ Vcb\, \ Vub/Vcb\, £k and 
AMa^s- On the other hand, the enhanced function C and its large complex phase can, in 
principle, affect the analysis of the UT through the double Z°-penguin contributions to 
Ek and AMs,d- However, we have checked that these effects are at the few-percent level. 
In particular, the phases of the B^ s^-^d s mixing amplitudes are essentially unchanged. 

Next the decay amplitude of the "golden" mode B^ — > J/ ipK^ is only insignificantly, 
i.e. at the level of the current experimental and theoretical uncertainties, affected by the 
NP in the EW penguin sector in our scenario. Consequently, as discussed in detail in 
Subsection 6.2, we may still convert - to a good approximation - the mixing-induced CP 
asymmetry AQp{Bd — > J/ipKs) into the angle P of the UT, as in the SM. In particular, 
a "universal unitarity triangle" [40] can also be constructed in our NP scenario by using 
only \Vcb\, AMd/ AMs and sin 2/3. This simplifies our phenomenological analysis signifi- 
cantly, since the determination of the CKM parameters can be separated from the study 
of the hadronic and NP effects in non-leptonic B decays and rare decays. However, as 
advertised in Step 2 of our approach, one can also determine 7 by complementing the 
B — ^ 7171 system with those B nK modes that are insensitive to the NP effects dis- 
cussed here, i.e. to EW penguin contributions, and subsequently compare the result with 
the one obtained by using the conventional UT fits. Following these lines and using the 
experimental information on the side Rb of the UT, we may also determine the remain- 
ing two angles a and (3 in an alternative manner, and eventually arrive at a remarkably 
consistent overall picture. 

Most interesting for us are the functions X, Y and Z, which are given by 

X{v)=C{v)+B''''{v), Y{v) = C{v) + Bf'!'{v), Z{v) = C{v) + ^D{v), (2.5) 

with C(v), D{v), B'^^{v) and B^^P'{v) resulting from Z'^-penguin diagrams, 7-penguin 
diagrams and AF = 1 box diagrams with uu and /i/i, respectively. Moreover, explicit 
calculations indicate that in the case of AF = 1 box diagrams with uu and dd in the 
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final state we have, to an exceUent approximation: 



(2.6) 



so that the functions X and Y in (2.5) are valid for non-leptonic decays as well. Similarly, 
it is found that the NP contributions to these AF = 1 box diagrams and to D are 
rather small, so that we may use for them, in the following analysis, their SM values 
{rrit — 167 GeV), which are given by 



As a consequence of these approximations, the dominant NP effects come from the Z^- 
penguin function C{v). In the standard MFV scenarios, also the one considered in [16], 
the function C{v) is real- valued. Here it contains a weak phase 9c- Consequently, the 
functions X, Y and Z are now given as follows: 



X(v) = |C(^;)|e*^^ + 0.73, Y{v) = \C{v)\e'^^ + Z{v) = \C{v)\e'^^ -Q.12. (2.8) 



While our analysis does not rely on a particular model with the properties specified 
above, concrete models with enhanced CP-violating Z°-mediatcd FCNC couplings gen- 
erated either at the one-loop level or even at the tree level have been discussed in the 
literature. They are reviewed in [31]-[34], in particular in the last of these papers; see 
also [43] . Also models with Z'-mediated FCNCs could be put in this class, provided their 
contributions can effectively be absorbed in the function C{v). For a recent analysis, 
see [44]. 

2.3 Relation Between Rare Decays and B tzK Modes 

As already discussed in [16], the connection between the rare decays and the B — > nK 
system is established by relating the function C to the EW penguin parameter q by 
means of a renormalization-group analysis. In the case of a complex C, the relation 
given in [16] is generalized as follows: 



where q and ^ characterize the EW penguin sector and enter the parametrization of the 
B 7iK decays given in Section 4. The numerical coefficients in (2.9) correspond to 
as{Mz) = 0.119 and depend very weakly on its value. 

In order to define q and 0, we recall formula (43) of [13] for the ratio of the EW 
penguin amplitude P^-y^ to the tree contributions in S — > irK decays (see also Section 4) , 
yielding 



are the Wilson coefficients of the {V — A)<^{V — A) EW penguin operators Qg and Qio, 



S""(i;) = B'^^'iv) = 0.73, B'^'^iv) = B^^'^^v) = 0.18, D{v) = -0.48. 



(2.7) 



|C(i;)|e^^c = 2.35 ge^"^- 0.82, q^q 



Vub/Vcb 

0.086 



(2.9) 
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respectively, which enter the effective Hamiltonian for AB = 1 non-leptonic decays [35], 
and the coefficients C[ 2{l^b) are given as follows 



C'M = Ci(//b) + \c,{ii,), C'M = C^ifx,) + ^CioM, (2-11) 

with Ci^2 being the Wilson coefficients of the current-current operators (5i,2- As Cg^io = 
0{a) and Ci^2 = we will approximate C{ 2 by Ci^2- Moreover, we will set = 1 

in front of Ciq. This approximation is justified by the fact that Cio{Mw) = and 
consequently Cio(/i6) is significantly smaller than C^{^b). Next, in the NP scenario 
outlined above, we may write 

C^il^h) + Cio(/x;,) = h + gC{v)e''^ = -^Ewe^^ (2.12) 

where h and g arc calculable by means of the two-loop renormalization-group formula 
given in [47]. The minus sign in (2.12) causes the parameter Agw to be positive in the 
SM. Writing then ^ = (^e*"^, r being a strong phase, we obtain 

CM + C2(/i6)ee- = Bcce-'^, (2.13) 

where it was convenient to introduce a strong phase ui [13], which vanishes in the SU{3) 
limit [45, 48]. Consequently, (2.10) yields the following simple expression: 



P' 



T' + C 



AS'^^-Sr^^c) ^ ,e^0e^^, q = ^tttt^^: (2.14) 

2 A|\/ub/V/c6| -Dec 



which we will use for the parametrization of the B — > irK decays in Section 4. 

As far as the rare decays considered in Section 5 are concerned, they can be directly 
expressed in terms of the X, Y and Z functions introduced in (2.8). Consequently, in 
order to discuss the correlation between the B — > ttK decays and the rare decays, it is 
useful to express them in terms of {q, (f)) by inserting (2.9) into (2.8). The corresponding 
expressions are given in Section 5. 



2.4 CKM Parameters 



Concerning the CKM parameters, we will use the Wolfenstein parametrization [49], gen- 
erahzed to include higher orders in A = ]T4s| [50]. Writing then 



Vtd = ARtX'e-"^, Vts = -\Vu\e- 
with tan (3s ~ —X^fj, we have 

Xt = V*Vtci = -fA^X^Rte-^'^e''^' with f = 

and 



Vts 



cb 



0.98, 



) = V*Vu = AX'Rte-'^, aJ^^ ^ VW^s = -fAA^e^^^ , 



(2.15) 

(2.16) 
(2.17) 
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where Rt is one of the sides in the UT, and A = iKbl/A^ is the usual Wolfenstein 
parameter. Moreover, we have the following useful relations for the two UT sides Rt and 
Rb and the ratio | Kb/Kb I [3]: 



Rf 



sm 7 



Rh 

















Vcb 


; 


Vcb 


-{ 



sin(/3 + 7) 
For our numerical analyses, we will use 

A = 0.2240 ± 0.0036, A = 0.83 ± 0.02, 



A 



sin P 



1-Xy2j sin(/5 + 7) 



. (2.18) 



V. 



ub 



0.086 ±0.008, = 0.37 ± 0.04 



/3 = (23.5 ± 2.0)°, (5s = -l°. 



(2.19) 
(2.20) 



2.5 Summary 

Before turning to the S — > tttt system, let us summarize the main results of this section: 

• The NP scenario considered here involves two parameters, C{v) and 9c- 

• The relevant EW penguin parameters for the B — > ttK decays are g, and uj. 

• The parameter sets (C, Oc) and (g, 0) are related through (2.9), which allows us to 
investigate the correlations between rare decays and non-leptonic B decays. 

3 The B TTTT System 

In the literature, S — > tttt decays are usually considered in the context of the determi- 
nation of the UT angle a [51, 52]. Here we shall use these modes from a very different 
point of view, which is inspired by the analyses performed in [53]-[56]. Very recently, 
the decay B^ tt+tt" was also discussed in [57] in the context of bounds on 7 and the 
UT. However, these approaches differ significantly from the one presented in [12] and 
here. In particular, the channels involving neutral pions are not discussed at all in [57]. 



3.1 Basic Formulae 

In order to make our analysis more transparent, we shall first neglect the EW pen- 
guin contributions to the B ^ tttt decays, which play a very minor role [58, 59]; using 
the isospin flavour symmetry of strong interactions, we include these contributions in 
Subsection 3.6, following [13, 60]. We may then write the S — > tttt amplitudes as 



V2A{B+ - 
A{B', - 
V2A{B', 



TT TT 



-[T + C] = -[T + C] 

TT+TT") = -[f + P] 

.7rV) = -[C-P], 



(3.1) 
(3.2) 
(3.3) 
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with 



P 
f 

C 



X'ARbe'"'[C + {Vtu-S)]. 



(3.4) 
(3.5) 
(3.6) 



The parameters of the CKM matrix entering these formulae have been defined in Subsec- 
tion 2.4, whereas the Vq describe the strong amphtudes of QCD penguins with internal 
g-quark exchanges (g G {t, c, u}),^ including annihilation and exchange penguins; T and 
C are the strong amplitudes of colour- allowed and colour-suppressed tree-diagram-like 
topologies, respectively; and S denotes the strong amplitude of an exchange amplitude. 
In the usual notation employed in the literature (see, for instance, [58]), the colour- 
allowed and colour-suppressed tree-diagram-like amplitudes 



T 

C 



X^AR^e^'^T 

X^AR^e'^C 



(3.7) 
(3.8) 



appear; they differ from the T and C amplitudes through the {Vtu — ^) pieces, which 
may actually play an important role, as was emphasized in [64]. However, we observe 
that these terms cancel in 

f + C = T + C, (3.9) 
and that the amplitudes in (3.1)-(3.3) satisfy the well-known isospin relation [51] 



V2A{B+ 7r+7r°) = A(S° ^ tt+tt") + V2A(S° 
It is convenient to rewrite these amplitudes as 



77 77 ). 



V2A{B+ - 
MB', - 
y/2A{B', 



■ 7r+7r°) = 
7r"'"7r~) = 



7t\') 



\P\e'' 



with 



C 
T 



-— e 
T 
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i{5p-5j,) 



1 + ^e*7e'(A-(^) 
d 



r-{Vtu-s) 

Vtc 



1 

Rb 



r-{Vtu-s) 



(3.10) 

(3.11) 
(3.12) 

(3.13) 

(3.14) 
(3.15) 



where the latter parameter was already introduced in [53], and 5f and 5p denote the 
CP-conserving strong phases of the amplitudes T and P, respectively. If we now consider 

^Strictly speaking in the case of Vt this standard terminology could be misleading, as the top-quark 
effects in amplitudes for B decays are only present in the relevant Wilson coefficients - also absorbed 
in Vq - and not in hadronic matrix elements. For relations between the language used here to that of 
the operator product expansion, see [61]-[63]. 
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the corresponding CP-averaged decay-amplitude squares, which are generically defined 
through 

{\A\')^\[\A{B^f)\'+\A{B^f)\'], 



we obtain 



{\A{B 



± 



[l + 2x cos A x^] 



(|^(Sd^7r+7r-)|^) = |f|2 [l-2dcos^cos7 + cZ^] 



TT TT 



1 + 2 (|^) cos( A - ^) cos 7 + 



(3.16) 

(3.17) 
(3.18) 

(3.19) 



These quantities provide two independent ratios of CP-averaged branching ratios, which 
we may choose as those introduced in (1.1) and (1.2). Using (3.17)-(3.19), we obtain 



1 2a; COS A 
1 — 2d COS 9 cos ^ + (P 



-'^00 



+ 2dx cos (A — 6) cos ^ + x"^ 



(3.20) 



(3.21) 



1 — 2d cos 9 cos 7 + (i^ 

In addition to and _Rqq , also the time-dependent CP asymmetries of the decays 
Ed — > 7r+7r~ and 3^ tt^tt" provide valuable information. Thanks to the efforts at the 
B factories, experimental results for the former channel are already available, which is a 
decay into a CP-even final state, exhibiting the following rate asymmetry [2]: 



r(i?°(t)^7r+7r-)-r(B»(t)^7r+7r-) 
= ^^^(5d ^ TT+TT") cos(AM,t) + (^rf ^ TT+TT") sin(AM,t). 



(3.22) 



Here, AM^ > is the mass difference of the mass eigenstates of the B^-mBson. system, 
while the CP-violating observables 



1 - le^? 



and (Ed ^ TT+TT" 



2Ime^? 

^ TT+TT 

1 + leit- 



(3.23) 



describe "direct" and "mixing-induced" CP violation, respectively, and are governed by 



id) 



—e 



-17 



-de' 



(3.24) 



The quantity (f)^, which equals 2/3 in the SM, is the CP-violating weak B^-B^ mixing 
phase. Consequently, we arrive at the following expressions: 



Acp{Bd 



TT TT 



2d sin 9 sin 7 



TT^TT 



1 — 2dcos^cos7 -I- 

sin(0rf -I- 27) — 2d cos 9 sin(0d + ^) + d'^ sin (pa 
1 — 2d cos ^ cos ^ + 



(3.25) 
(3.26) 
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The current experimental status of these obscrvables is given by 

idir.o -0.19 ±0.19 ±0.05 (BaBar [65]) 

^CPl^d^TT TT J - |_o.77±o.27±0.08 (Belle [66]) ^'^'^'^ 

/imix.R _ + -^ _ / +0.40 ± 0.22 ± 0.03 (BaBar [65]) 
Acp{Bd^n TT +1.23 ± 0.41^°:°^ (Belle [66]). ^'^■^^> 

Unfortunately, the BaBar and Belle results arc not fully consistent with each other, 
although both experiments point towards the same signs, and the last BaBar update of 
AQp{Bd TT'^TT") has moved towards Belle. In [11], the following averages are given: 

Acp{Bd Tr+TT") = -0.38 ± 0.16 (3.29) 
Acp(Bd TT+TT") = +0.58 ± 0.20. (3.30) 

As was pointed out in [55, 56], the CP asymmetries in (3.29) and (3.30) point towards 

7 ~ 60° for (f)d ~ 47°, in accordance with the SM. Before discussing the determination of 
7 in more detail in Subsection 3.4,^ it is useful to first have a closer look at the hadronic 
5 — > TTTT parameters. 

3.2 Determination of the Hadronic Parameters 

In order to explore the hadronic S — > tttt parameters, we assume that 

7 = (65 ±7)°, 0rf = 2/3=(47±4)°, (3.31) 

as in the SM [4]. If we then look at (3.25) and (3.26), we observe that each of these 
observables allows us to determine d as a function of the strong phase 9. In the case of 
the direct CP asymmetry Acp{Bd — > tt'^tt"), we obtain 



d 



Aqp {Bfi ^ tt~^7t ) COS 9 cos 7 — sin ^ sin 7 



A^I.{Ba^n+n-) 
±\l [A'^l,{Bd TT+TT") COS ^ cos 7 - sin 9 sin 7] ^ - A'^l>{Bd tt+tt")^ 
whereas its mixing-induced counterpart A^p{Bd tt+tt^) inipUes 



(3.32) 



d^k± - /, (3.33) 

with 



k 



sin(0rf + 7) - A'^^jBd TT+TT-) cos 7 

Sin0rf-^g^^-(fid^ TT+TT-) 



cos^ (3.34) 



; ^ sin(0, + 27)-^bT(i?^^TT+TT-) 

sin0rf-^-^?(5,^TT+TT-) ■ ^ ' ' 

In Fig. 1, we show the corresponding contours, and observe that we obtain a twofold 
solution for [d,9). It should be emphasized that these contours and the corresponding 
determination of (d, 9) for a given value of 7 are theoretically clean. 

^See also Appendix B, where a very recent update by the Belle collaboration [67] is discussed. 
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50 100 „ 150 200 250 



Figure 1: The contours in the 6-d plane for 7 = 65° and (pd = 47°. The sohd hues 
correspond to the central value and la upper and lower ranges of Aqp^B^ tt+tt^) = 
-0.38 ± 0.16, the dashed lines represent A^'^{Ba n+n-) = +0.58 ± 0.20, and the 
dotted lines refer to = 7.17 ± 0.75, as discussed in Subsection 3.4. 

Let us now also consider the observables and -Rqo • complement them with 

the CP-violating observables AQp{Bd — >■ vr+vr") and Aqp^Bci — > tt+tt"), which allow us 
to determine d and 9 as we have just seen, we may extract the parameters x and A as 
well. To this end, it is convenient to introduce 

Rl'[_ = DRIZ_ (3.36) 
i?^o" = DR-,l (3.37) 

with 

D = 1 -2c/ cos ^ cos 7 + ^2. (3.38) 
Using then (3.20) and (3.21), we obtain 

X = - cos A ± - sin^ A (3.39) 

and 

x = -d cos 7 cos( A -e)± sJri^ - [1 - cos2 7cos2(A - 6)] d^, (3.40) 

respectively, allowing us to calculate two contours in the A-x plane. The intersections 
of these curves then allow us to extract x and A. Interestingly, the {d, 6) solution with 
d ~ 5 gives only complex solutions and can therefore by excluded. On the other hand, 
in Fig. 2, we show the contours in the A-x plane corresponding to the physical d ~ 0.5 
case, yielding a twofold solution for (x. A). 

Using (3.31) as an input, the expressions in (3.20), (3.21), (3.25) and (3.26) allow us 
to convert the data in (1.1), (1.2), (3.29) and (3.30) into the hadronic parameters {d,9) 
and (x. A). Following these lines, we obtain 

ci = 0.48l°l, ^ = + (138+1^)°, (3.41) 
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Figure 2: The contours in the A-x plane. Using d = 0.48 and 6 = 138° following from 
the central values of Fig. 1, we obtain the solid set of contours for = 2.12 ± 0.37 
and the dashed set of contours for Rqq = 0.83 ± 0.23. 

as well as the twofold solution 

a; = 1.22;°i?, A = -(71^^^)°, (3.42) 

x = imtril, A = + {52tlt)\ (3.43) 

where our treatment of errors is as described in Appendix C. At this stage, we cannot 
distinguish between the two solutions for xe*^. However, as we will show in Subsec- 
tion 4.4.3, the solution in (3.43) can be excluded through the B — > irK data, as it would 
correspond to large direct CP violation in tt^K^, which is ruled out by experiment. 

For the time being, we will hence focus on (3.42). Let us finally note that the determina- 
tion of the hadronic parameters given in (3.41)-(3.43) is essentially theoretically clean, 
and that the experimental picture will improve significantly in the future. 

3.3 Interpretation of the Hadronic Parameters 

The result of x = 0{1) obtained above, which implies that \C\ ~ |T|, is in stark contrast 
to the naive expectation of \C\ ~ 0.25 x \T\. At first sight, this feature seems to imply a 
complete breakdown of the concept of colour suppression in B ^ nn decays. However, 
one should realize that the usual arguments leading to \C\ ~ 0.25 x |T| neglect the con- 
tributions iVtu^^) in (3-5) and (3.6), which may be important [64]. In fact, our analysis 
suggests that these contributions cannot be neglected if we want to understand the pat- 
tern of the B ^ TCTT observables in a plausible manner. Because of the different signs in 
(3.5) and (3.6), we may explain the surprisingly small Bd — >■ vr+vr" branching ratio natu- 
rally through destructive interference between the T and {Vtu — £) amplitudes, whereas 
the puzzling large B^ tt'^tt^ branching ratio originates from constructive interference 
between the C and {Vtu — £) amplitudes. 
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On the other hand, we have 01^^^. = 180°. Consequently, as can be seen in (3.12), 
Ba — > 7r+7r~ would favour - in contrast to the SM expectation - 7 > 90° within the fac- 
torization approach, since BR{Bd — > tt'^tt') would then be reduced through destructive 
interference between trees and penguins. In contrast, we arrive at a picture with large 
non-factorizable contributions, exhibiting certain interference effects at the hadronic level 
that allow us to accommodate straightforwardly any currently observed feature of the 
S — > TTTT system in the SM. 

3.3.1 Insights into the Substructure of xe^^ 

The calculation of the hadronic parameter xe^^ from first principles is extremely chal- 
lenging and cannot be done in a reliable manner because of our limited knowledge of 
non-perturbative strong-interaction physics. However, in order to obtain deeper insights 
into the substructure of xe*^, it is instructive to write 

■KIT iA^^ I ^ iAf 



with 



are^^^'' ^ I (3.45) 



JAr — ^tU - ^ 



Consequently, we may convert xe*"^ into the complex quantity (e^^^ with the help of 

Ce'^ = ^ X -0 ■ 3.47 

^ l + 2xcosA + x2 ^ ^ 

In the special case of a'^'^e^'^^'' = xe^^, this expression reduces to Ce*^*^ = 0. However, a 
value of 

^^TTgiAr _ ^ X e-^^i° would appear as completely unrealistic. On the other 
hand, in view of the large non-factorizable effects exhibited by (3.41) and (3.42), we think 
that aY may well take values as large as 0.5, with a large strong phase ■ Fig- 3, 
we consider the central values of x and A in (3.42), and show the contours in the A^-^ 
plane corresponding to different values of d^'^ ^ where each point is parametrized by the 
value of /^Y e [0°,360°]. We observe that the "naive" value of a^^e^^a" ~ 0.25 would 
result in a rather large value of C ~ 0.65. On the other hand, a^e^'^^'' ~ 0.5 x e*^^°° gives 
a significantly smaller C, ~ 0.4, which may well originate from constructive interference 
between the Vt and Vu amplitudes, satisfying, for instance, \Vt/T\ ~ \Vu/T\ ~ 0.25. 
Moreover, for such values of d^e^^'^'^ , we have |1 + oYe^^^^'l ~ 1.25, in accordance 
with the "naive" expectation. Consequently, in the B"^ tt^tt^ channel, which is 
measured in agreement with the conventional theoretical expectations, not only the 
hadronic interference effects originating from the {Vtu ~ amplitude cancel, but also 
the large non-factorizable effects in aYe-^^^" would not manifest themselves. 

We shall come back to these considerations in Subsection 4.2.3, where the quantity 

CV^ . = M (3.48) 
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A. 



Figure 3: The contours in the A^-^ plane corresponding to the central values of xe*^ in 
(3.42) for various values of aj'^ between 0.2 and 0.5 and /^^ ^ [0°,360°]. The two dots 
refer to cases discussed in the text. 
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Figure 4: The contours in the A^;-^ plane corresponding to the central values of xe*^ in 
(3.42) for various values of between 0.2 and 0.5 and e [0°, 360°]. 
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will play an important role in the context of a hadronic parameter entering the charged 
B — > -kK decays. In Fig. 4, we show - in analogy to Fig. 3 - the contour plots in the 
A^-C plane for given values of and £ [0°,360°]. Interestingly, the analysis of 
Pc will also point towards a^e^^^'' ~ 0.5 x e'^^°°, thereby complementing the picture 
described above. 

Let us finally note that the hadronic parameters in (3.41) and (3.42) also imply 

= (0.27l°:13) xe^(-2i^3r, (3.49) 
yielding 

^^ = (0.10lH^)xe^(-^^^)°, (3.50) 

where we have used the value of Rb in (2.19). Consequently, the small numerical value 
in (3.50) would then require destructive interference between the Vt and Vc amplitudes, 
where Vc contains "charming" penguins [68] , in contrast to the constructive interference 
between Vt and Vu- Interesting insights into these subtle dynamical issues can be ob- 
tained with the help of the penguin-induced decay — > K^K^ [69], which can nicely 
be complemented with the Bg —>■ K^K^ channel [2]. 



p 


gn = _ 




1 


Vtc 


T + C 




1 + xe'^ 


Rb 


T + C 



3.3.2 Probing Penguin Annihilation and Exchange Topologies through the 
Bd K+K-, TT+TT- System 

The origin of the value of xe^^ in (3.42) could also be due to anomalously enhanced 
penguin annihilation and exchange topologies (which contribute to the B —>■ ttt: modes 
but not to their B — > ttK counterparts discussed in Section 4). However, the importance 
of these topologies can be probed through the decay B^^ — > K'^K~ , where the current 
experimental upper bound is given as follows [11]: 



BR{Bd K+K-) < 0.6 X 10"^ (90% C.L.). 



(3.51) 



The corresponding decay amplitude has the same structure as (3.12), with the impor- 
tant difference that now only penguin annihilation and exchange topologies contribute. 
Employing the SU (3) flavour symmetry of strong interactions, we may hence write 



A{B', ^ K+K-) = -X'AR, [S - {VA)tu] [e'^ + QvAe'^^^] , 



with 



1 

Rb 



(VA) 



tc 



(3.52) 



(3.53) 



£-{VA)tu. 

where the penguin annihilation amplitudes {VA)tu and {VA)tc are implicitly included 
in the B —>■ ttt: amplitudes Vtu and Vtc, respectively. Consequently, we obtain 



BR{Bd K+K-) 
BR(S± ^ 7r±7rO) 



S - {VA)tu 



T + C 



^1 2q'pj, cos d-pA COS 7 + qI,j^ ^ 0.2, 

(3.54) 
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which does not indicate any anomalous behaviour. In the future, this bound can be 
improved significantly. The decay — > K'^K~ is very accessible at LHCb, where it 
may be possible to reach the 10""^ level for its CP-averaged branching ratio [70]. At this 
experiment, we may also exploit the physics potential of the — > tt"*"??" channel, which 
is the [/-spin partner of the i?J] K^K^ modc,^ and has a decay amplitude of the same 
structure as the 5° — > K^K~ mode [53] (or the B'^ — > ti' transition discussed in 
Section 4). If we then use the [/-spin flavour symmetry and consider the CP-averaged 
branching ratios, we obtain 



VA 



BR(5rf ^ K+K- 



BRfS, 



TBI 



1 -|- 2qj,j\^ cos i^p^ cos 7 



2eg-pA cos -d-pA cos 7 + 



where we neglect tiny phase-space efi'ects, and 



e = 



1- A2 



= 0.053. 



(3.55) 



(3.56) 



The quantity H-pj^^ allows us to obtain constraints for the parameter g-pj,. Since the 
amplitude structure of the B^ — > K~^K^ , Bg — > 7i~^7i~ system is analogous to that of the 
Bfi — > TT+TT", Bs — > K'^K~ {Bd — > TT^K^) system, we may apply the bounds derived in 
[54], yielding 



1 — e^/Hp^ 



< QVA < 



1 + e^Hpjs, 



(3.57) 



1 - \jBpj^ 
> K^K~ strategy proposed in [53], a 



TT ' TT 



observables would allow a 



In complete analogy to the B^ n^7i~, Bg 
measurement of the CP- violating B^ K^K^, Bg 
determination of both 7 and the hadronic parameters gp^ and i^p^- However, the main 
interest in an analysis of the B^ K^K~ , Bg n'^n' system is obviously to obtain 
insights into the importance of penguin annihilation and exchange topologies. Using 
the corresponding information on and '&pa, the expression in (3.54) allows us to 
determine \{S - {VA)tu)/{'T + C)|, while 



BR{Bg TT+TT-)' 

BR(S± ^ 7r±7rO) 



1 

Rb 



{VA)tc 



T + C 



(3.58) 



It will be very interesting to confront these considerations with experimental data for 
the Bg — > 7r"^7r~ channel. First constraints may soon be available from run II of the 
Tevatron, which should be considerably improved at LHCb. 



3.4 Determination of 7 



Let us now come back to Fig. 1. There we have also included a contour, which is related 
to the following quantity [54] : 



H 



c 


2 


'c 





Mb, ^Mk/Mb^,Mk/MbJ r^o 
Mb^ ^MJMb.^MJMb,) Tbo 



BR{Bd TT+TT") 

BR{Bg K+K-) 



'The [/-spin flavour symmetry of strong interactions connects the down and strange quarks in the 
same manner as the ordinary SU{2) isospin flavour symmetry connects the down and up quarks. 
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U spin 1 — 2d cos 9 cos 'f + (P 
+ 2ed cos 6 cos 7 + c?^ ' 

where \C'/C\ is a t/-spin- breaking parameter, and 



(3.59) 



$(x, y) = y/[l-{x + yy] [l-ix-yf] (3.60) 

denotes the usual two-body phase-space function. The introduction and the use of H, 
which has the same structure as (3.55), is inspired by a variant of the strategy presented 
in [53] , allowing us to determine 7 and the hadronic parameters {d,9) by relating the 
CP- violating — > 7r"'"7r~ observables to their Bg — ^ K^K^ counterparts through the 
[/-spin flavour symmetry of strong interactions (see also Subsection 4.6.1). This strategy 
is very promising for B-decay experiments at hadron colliders. First important steps 
are already expected at run II of the Tevatron [71], whereas it can be fully exploited 
in the era of the LHC [72], in particular at LHCb, where experimental uncertainties 
of 7 at the few-degree level may be achieved [73]. Unfortunately, the physics potential 
of Bg K^K~ cannot yet be exploited. However, if we assume that the penguin 
annihilation and exchange topologies discussed in Subsection 3.3.2 play a minor role 
and employ SU (3) flavour-symmetry arguments, we may replace B^ — > K^K^ through 
Bd — ^ TT^X^, and arrive straightforwardly at 



2 



BR{Bd TT+vr-) 



7.17 ±0.75, (3.61) 



where the ratio fx/ fw = 160/131 describes factorizable 5'C/(3)-breaking corrections [54]. 
In fact, the utility of B^ — ^ n'^K^ transitions to address the "penguin problem" in 
Bfi ^ vr+vr" was already pointed out about ten years ago [74]. If we now use the 
expression for H in terms of the hadronic parameters given in (3.59), we obtain [54]: 

d——b cos 9 cos 7 ± ^/a -I- (6 cos 9 cos 7)^, (3.62) 

with 

l-e'H 1 + eH 

a = ^-^, b=^^. (3.63) 

which allows us to calculate d as a function of 9 for 7 and H as given in (3.31) and 
(3.61), respectively. The resulting contour is shown in Fig. 1 and fits perfectly into the 
picture of the theoretically clean 9-d contours following from the CP-violating B^ 
7r+7r~ observables. The fact that the H contour described by (3.62) goes right through 
the intersection of (3.32) and (3.33) is very remarkable and indicates that the penguin 
annihilation and exchange amplitudes as well as non-factorizable SU (3)-breaking effects 
actually play a minor role. Otherwise, the corresponding hadronic uncertainties would 
have to conspire in a very contrived way, which does not look plausible to us. 

If we do not fix 7 through the UT fits, but use instead H, as discussed in [53]-[56], 
we may determine this angle, yielding the following twofold solution: 

7 = (40.5l|l)° V (64.711^)°. (3.64) 
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The errors are obtained with the help of a standard analysis from the errors of 
AQp{Ba — > 7r+7r~), AQp{Ba tt+tt") and H by fitting to the hadronic parameters and 
7. In addition, there is an error of ±2.5° on 7 from the uncertainty in (pd- A detailed 
analysis of the whole B — > ttK system, as discussed in Subsection 4.3, allows us to 
lift this degeneracy, leaving us essentially with the solution around 65°, which agrees 
perfectly with the range for 7 given in (3.31). 

An interesting alternative to determine 7 is provided by the direct CP asymmetry 

^ BR{B^ TI-K+) + BR(50 ^ n+K-) ^ ' 

If we employ again the SU (3) flavour symmetry and the plausible dynamical arguments 
specified above, we may express this observable in terms of the hadronic parameters 
introduced in (3.15) as follows: 

_ ^ . , 2e(isin^sin7 
-^°^'^^--"^) = e^ + 2edccs9ccsT + d^ - (3.66) 

In contrast to expression (3.61) for H, (3.66) does not involve a {fx/ IttY term and is 
hence affected to a smaller extent by S'?7(3)-breaking corrections, i.e. is more favourable 
from a theoretical point of view. On the other hand, as can be seen in (3.65), a non- 
vanishing value of AQp{B(i tt^K^) has not yet been established by the B factories, 
although the experimental evidence for this CP asymmetry is steadily increasing. 
In this context, it is important to note that we have the following relation [53, 54]: 



A'^^iBd ^ TT+TT-) 



fK 

U 



0.25±0.13 0.38±0.04 



where we have also indicated the current experimental results. Within the uncertainties, 
this relation is satisfied by the data. In particular, also the pattern of the signs of 
the direct CP asymmetries is in accordance with (3.67). Should the central values be 
correct, they would correspond to non-factorizable 5'[/(3)-breaking corrections at the 
B — > TTTT, ttX amplitude level of 0(20%). On the other hand, the current central value 
for AQp{Bd —>■ TT^K^) may also be on the lower side, which would be our preferred 
scenario. It will be exciting to observe how the data evolve in the future. 

In our analysis, we have assumed that 4>d ~ 47°, in agreement with the SM. However, 
as discussed in [55, 56], it is interesting to consider also the second, unconventional 
solution of (f)d ~ 133°. If (pd, 7, d, 0, x, A are solutions of (3.20), (3.21), (3.25), (3.26), 
(3.59), (3.66), then 

TT — (pd, — li d, TT — 9, X, —A (3.68) 

are solutions as well. Consequently, (3.68) allows us to go easily from the (pd ~ 47° to the 
(pd ~ 133° case. Interestingly, for the value of 6 in (3.41), wc obtain cos^ ~ —0.7 < 0, 
exhibiting the same sign as within factorization, where ^If^c^ = 180°. On the other hand, 
the value of 9 corresponding to (pd ~ 133° yields cos^ ~ +0.7 > 0, i.e. the opposite 
sign, thereby disfavouring the 0^ ~ 133° solution. Moreover, also the experimental 
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Figure 5: The dependence of AQp{Bd tt^tt") on 7, as described in the text. 

sign information about CP violation in Bd D^*^^!^^ modes, obtained recently by the 
BaBar collaboration, points towards (pd ~ 47° [75]. Unfortunately, the Belle collaboration 
favours the opposite sign [11]. 

We shall return to A'^p{Bd 'k^K'^) in Subsection 4.4.1, and to the determination 
of 7 in Subsection 4.3, where we shall also address the other two angles, a and /3, of the 
UT. For the following numerical analyses, we will continue to use 7 as given in (3.31). 

3.5 Prediction of CP Violation in Bd tt^tt^ 

Let us consider, as the last element of our B ^ titt analysis, the CP-violating observables 
of the decay Bd n^ir^, which is - in analogy to Bd — >■ vr+vr" - another transition into 
a CP-even final state. If we apply the standard formalism as discussed for Bd —>■ tt^it~ 
in (3.22)-(3.24), we obtain 

/Idir/R 2dxsin(g- A)sin7 

AM - vr ) = ^. + 2da:cos(^-A)cos7 + x^ ^^"^^^ 

■ (P sin ^d + '^dx cos{9 - A) sin(0d + 7) + sin((^^ + 27) 

Acp [Bd ^ n n ) = — — . (3. ([)) 

^ c?2 + 2c?xcos(6'- A)cos7 + x2 ^ ^ 

Complementing these formulae with the analysis outlined above, we may predict these 
CP- violating observables, with the following SM result: 

A'S^AB, - = -OAlt'of, (3.71) 

A^^iBd ^ rr%%^ = -0.551°:^, (3.72) 

where the errors come from the procedure described in Appendix C. In the future, 
when the errors of our input quantities will decrease, more accurate predictions for these 
quantities will be possible. 

It is interesting to note that a measurement of one of the CP asymmetries of the 
Bd tt^tt" mode would allow an essentially clean determination of 7. In Figs. 5 and 
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Figure 6: The dependence of AQp{Bd tt^tt") on 7, as described in the text. 



6, we illustrate the dependences of Aqp^B^ — > 7r°7r°) and AQp{Bd — >• tt^tt") on this 
angle, respectively, where we consider the central experimental values of K^_, Rqq, 
AQp{Bd vr+vr^) and Aqp^B^ vr+vr^) for simplicity.^ We observe that in particular 
AQp{Bd 71^71^) - being more sensitive to 7 than AQp{Bd ti^ti^) - would allow a 
useful determination of this angle, which would also be unambiguous in the range of 7 
considered. 



3.6 EW Penguin Contributions 

Let us finally come back to the EW penguin contributions to the B — > vrvr modes, which 
we have neglected so far. Using the results of [13, 60], we may take them into account 
with the help of the isospin symmetry of strong interactions. If we consider the SM and 
absorb the colour-suppressed EW penguins in the amplitude P, (3.1)-(3.3) are modified 
as follows: 



V2A{B^ 



V2A{B', 



TT TT 



TT TT 



,«7 



iPle" 



X 



1 + -e'^e' 
a 



+ q 



1 + xe 



(3.73) 
(3.74) 

(3.75) 



where 



q = 



Pt 



EW 



T + C 



1.3 X 10^2 X 



td 



V,, 



ub 



1.3 X 10^2 X 





sin 7 




sin [3 



3 X 10~l (3.76) 



The impact of the EW penguin contributions on the determination of the hadronic 
parameters {d,6) and (x. A) discussed above is negligible. For q as in (3.76), x changes 

^We suppress also certain unphysical branches in these figures, which are related to unphysical 
solutions for the hadronic parameters d, 0, x and A. 
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by about 1% and A by 1°, whereas d and 6 are not affected by q. Even if an arbitrary 
(NP) weak phase for q is admitted, x changes by at most ±4% and A by ±2°. 

In contrast to the S — > tttt decays, there are modes in the B nK system that 
are significantly affected by EW penguins [58, 59, 76], since the CKM structure of these 
channels is very different from the B tttt case. Despite complications for the explo- 
ration of the B — > ttK system, this feature offers a nice avenue for NP to manifest itself 
in the S-factory data. 

3.7 Summary 

Before having a closer look at the B — > ttK system in the next section, let us summarize 
the main results of our B ^ tttt analysis: 

• The starting point is a general parametrization of the relevant decay amplitudes. 
Assuming that 7 and (pa = 2/3 agree with the SM expectations, we may extract the 
hadronic parameters from the S — > tttt data in an essentially clean manner. 

• This analysis results in values of the hadronic parameters that show large non- 
factorizable contributions. In particular, we have seen that the amazingly large 
Bfi — > tt'^tt^ branching ratio and the surprisingly small B^ — > tt'^tt' branching ratio 
can be conveniently accommodated in the SM through constructive and destructive 
hadronic interference effects, respectively. The B^ vr^vr^ channel, which is 
measured in accordance with the conventional theoretical estimates, is not affected 
by this mechanism. Moreover, we expect that also the B ^ nn colour-suppression 
factor is sizeably affected by non-factorizable contributions, where, for example, 
Q^TTgiAj'^ ~ 0.5 X e~*™° would fit nicely into our picture. 

• The Bd — > K'^K~, Bg — > tt+tt" system allows us to probe the importance of 
penguin annihilation and exchange topologies, which are expected to play a minor 
role. This strategy is particularly promising for the era of the LHC The current 
5-factory bounds on B^, K^K~ do not indicate any anomalous behaviour. 

• Assuming that penguin annihilation and exchange topologies are negligible, we 
may complement the B — > tttt data in a variety of ways with the experimental 
information provided by B^ — * ti'^K^ modes, which are only marginally affected 
by EW penguins. Following these lines, we may obtain insights into SU (3)-breaking 
effects and may also determine 7. Interestingly, we obtain results in good agreement 
with the well-known UT fits, and arrive at a remarkably consistent overall picture. 

• Finally, having all relevant parameters at hand, we may predict the CP-violating 
observables of the decay B^ — > tt^^tt^, with the promising perspective of having 
large direct and mixing-induced CP violation in this channel. A measurement of 
one of the corresponding CP asymmetries would allow a clean determination of 7. 
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4 The B ttK System 

Decays of the kind B — >■ ttK have received a lot of attention in the hteraturc over the 
last couple of years (for reviews, see [2, 77]), mainly in the context of the determination 
of the UT angle 7. However, as was emphasized in [13, 59, 78], these modes offer also 
valuable insights into the world of EW penguins. Interestingly, puzzling features of the 
current S-factory data for the B — > nK modes may point towards NP effects in the EW 
penguin sector, as we will show in this section. A crucial ingredient of this analysis is 
given by the hadronic B tttt parameters determined in Section 3, which emerged from 
our resolution of the "5 — > tttt puzzle" within the SM. 



4.1 Basic Formulae 

If we employ the isospin flavour symmetry of strong interactions, we may decompose the 
B — > ttK amplitudes in the following manner: 

(4.1) 
(4.2) 
(4.3) 
(4.4) 



A{B+^ 


7r+X°) 


= -P' [I + pce'^-^e^^] 




V2A{B+ ^ 


tt'K+) 


= P'[l + p^e'^^^e'^ - (e*^ - 






TT-K+) 


= P'[l- re'^e'^] 




V2A{B', - 


> 7r°X°) 


= -P' [l+p^e'^-e'^ -qe' 





Here the CP-conserving strong amplitude 



P'^(l-^) AX\V[-V':) 



(4.5) 



is the B — > ttK counterpart of (3.4),^ describing the difference of the QCD penguins 
with internal top- and charm-quark exchanges. 



1^ 



v[-K-A' 

VI - V' 



(4.6) 



where 7^4 is the strong amplitude of QCD penguins with internal up-quark exchanges 
contributing to the charged B — > ttK decays and A' denotes an annihilation amplitude. 



X^Rb 
1- A2 



T + C 



(4.7) 



where T' and C are the colour-allowed and colour-suppressed tree- diagram-like topolo- 
gies corresponding to their B ^ tttt counterparts T and C in (3.5) and (3.6), respectively. 



re 



''The primes remind us that we are deahng with b ^ s transitions. 



X^Rb 
1^ 



(4.8) 
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where V!^ is the strong amphtude of QCD penguins with internal up-quark exchanges 
contributing to the neutral B ttK decays, 



X^Rb 



(4.9) 



and the EW penguin parameter qe^'^e^ was introduced in (2.14). In (4.1)-(4.4), we 
have neglected certain contributions from colour-suppressed EW penguins, which are 
expected to have a minor impact on our analysis. Since a detailed discussion of these 
topologies is rather technical, we have left it to Appendix D, where also the correspond- 
ing generalizations of (4.1)-(4.4) can be found. As we will see in Subsection 4.4, the 
current S-factory data for those B nK modes where EW penguins may only con- 
tribute in colour-suppressed form do not indicate any anomalous behaviour, i.e. do not 
point towards an unexpected enhancement of the colour-suppressed EW penguins. The 
amplitudes in (4.1)-(4.4) satisfy the following well-known isospin relation [79]: 

A{B+ ^ 7r+X°) + V2A{B+ ^ 7r°i^+) = A(5° ^ 7r-K+) + V2A(S° ^ 7r°X°) 

= - [e'^ - qe''>'e^] \T' + C'\e'^T'+^' , (4.10) 

as well as [59] 

Let us also note that the hadronic parameters pce^^", rcc'^^", re^^ and pnC*^" introduced 
above do not involve any EW penguin contributions, in contrast to the quantities appear- 
ing in the parametrization proposed in [13]. This feature is important for the following 
considerations. 



4.2 Determination of the Hadronic Parameters 
4.2.1 pce^^-^ 

As can be seen in (4.6), the parameter PcC*^" entering the charged B ttK decays is 
expected to be tiny because of X^Rb ~ 0.02. On the other hand, a sizeable value of Pc 
would be indicated by a significant direct CP asymmetry 

BR{B+ 7r+X°) - BR{B- n-R^) 



BR{B+ n+K^) + BR{B- n-R^) 



()^ sill 9,. sill 



1 -I- 2pc cos 9c cos 7 -I- p; 



and an enhancement of the CP-averaged -B^ — > K^K branching ratio [62, 
the current S- factory results for these quantities [11], 



BR{B^ 



^ TT^K) = -0.02 ± 0.06 
K^K) < 2.4 X 10-^ (90% C.L.), 



(4.11) 
However, 

(4.12) 
(4.13) 
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do not indicate any anomalous behaviour. In particular, if we employ the f/-spin flavour 
symmetry of strong interactions and introduce (for a detailed discussion, see [2]) 



K 



,p2 



BR(B± ^ K^K) 



1 + 2pc cos 6c cos 7 + p, 



— 2epc cos 9c cos 7 + p; 



,2' 
c 



(4.14) 



we obtain the following allowed range for pc, which has the same structure as (3.57): 



l-e^fK 



<Pc< 



1 + 



II - VK 



(4.15) 



Using now Rsu{3) — 0.7, which describes factorizable [/-spin-breaking corrections [2], 
and BR(S± ^ K^K)/BR{B^ ti^K) < 0.1, which follows from (4.13) and the 
experimental result BR(S± n^K) = (21.8 ± 1.4) x 10"^ [11], we arrive at 



Pc < 0.1. 



(4.16) 



In the future, this bound can be improved signiflcantly. We shall neglect the parameter 
Pc in the following discussion. Should K^K decays soon be observed at the B 

factories, thereby indicating a value of pc at the 0.1 level, this parameter could be taken 
into account by following the strategies discussed in [13, 46, 62, 80]. In this context, it 
should be emphasized that only the analysis of the charged B nK system may be 
affected by Pc, whereas this quantity does not affect the neutral B ttK decays. 



4.2.2 Vce"-^", re^" and pnC***" 

Let us now turn to the other hadronic parameters Vce^^", re^^ and pnC*^" appearing in 
(4.1)-(4.4). If we look at (4.7), (4.8) and (4.9), we may derive the following relations: 



" — re + PnC ° 

jen 



C' + VL 



tu 



(4.17) 

(4.18) 

Consequently, Tct^^" and re"^ differ through the quantity PnC*^", which is proportional 
to x'e^^ . Let us next assume that the penguin annihilation and exchange topologies 
discussed in Subsection 3.3.2 play a minor role. Using then the SU{3) flavour symmetry 
of strong interactions, we may relate the hadronic B nK parameter x'e^^ to its 
B nil counterpart xe"^^ introduced in (3.14) through the simple relation 



with 



xe 



iA 



(4.20) 



As far as the parameter re** is concerned, we have 



re 



d 



(4.21) 
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with e defined in (3.56). Consequently, (4.17), (4.18), (4.20) and (4.21) allow us to 
determine re**^, pnC*^" and rcC'^" from the S — > tttt analysis performed in Section 3. 
Following these lines, the numerical values given in (3.41) and (3.42) imply 

r = 0.1lt'oZ S = +{42tllr (4.22) 
Pn^O.lSt'oZ e,^-{29tllr (4.23) 

rc = 0.20lH?, 5c = +(211)°, (4-24) 

where our treatment of errors is described in Appendix C. Interestingly, the values in 
(4.24) mimic the picture of QCD factorization [9, 81], whereas (4.22) and (4.23) differ 

strongly from the corresponding predictions. 

The parameter Tc in (4.24), which follows directly from the simple expression 

Tc = ^ ^[1 - 2d cos e cos 7 + d^] Rl'L , (4.25) 
can be determined alternatively with the help of the following well-known relation [82] : 



V2 





fK 


/ BR(B± ^ 7r±7rO) 


Vud 


fn\ 


/ BR(5± ^ TT^i^O) 



0.196 ±0.016, (4.26) 



which relies on the SU{3) flavour symmetry and the neglect of the pc term in (4.1). We 
consider the agreement between the numerical values in (4.24) and (4.26) as very re- 
markable. In particular, (4.24) does no^ rely on any assumption about the Pc parameter, 
thereby supporting its neglect, in addition to the arguments given above. Whereas the 
impact of pc on (4.26) is maximal for ~ 0° V 180°, as this strong phase enters there 
through cos^c, (4.11) is governed by sin^c, and is hence affected most for 9c ~ ±90°. 
Consequently, we obtain complementary information on pc, suggesting that this parame- 
ter is indeed negligible. Interestingly, also enhanced colour-suppressed EW penguins (sec 
Appendix D) could affect the determination of Tc through (4.26), thereby leading to a 
possible discrepancy with the value following from (4.25). The agreement between these 
two determinations of Tc does therefore also not point towards an anomalous behaviour 
of the colour-suppressed EW penguins. Moreover, it suggests moderate non-factorizable 
S'f/(3)-brcaking corrections.^ A comprehensive study of the SU (3)-brcaking corrections 
to the relation in (4.26) was recently performed with the help of QCD sum rules [83], 
with the result that the dominant effects are actually described by the "factorizable" 
fx/ fir factor, which is in accordance with the picture following from our phenomeno- 
logical analysis. In this context, it is also interesting to note that (4.21) is not affected 
by any 5'[/(3)-breaking corrections within the factorization approximation, as the cor- 
responding decay constants and form factors cancel [53]. On the other hand, (4.20) is 
affected by factorizable SU (3)-breaking corrections in the following manner: 



xe'^, (4.27) 



^Otherwise, these effects fiave to conspire in a very contrived way with the PcC**" parameter and the 
colour-suppressed EW penguins, which does not seem plausible to us. 
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where the Bauer-Stech-Wirbel model [84] points to a small deviation from 1 of the cor- 
rection factor, which we shall neglect in the following discussion. These SU (3)-breaking 
effects are neglected in (4.26) as well. 



i0c 



4.2.3 Further Theoretical Aspects of p^e 

From the theoretical point of view, it is instructive to complement the considerations of 
Subsection 3.3.1 with our B nK results, and to explore the implications for PcS^^"- If 
we employ (3.46) and the relations given above, we obtain 



Pee' 




A' 



1-Ce 



re 



A' 



(4.28) 



where we have used the parameter (e^^( introduced in (3.48) and have taken the nu- 
merical values in (4.22) into account. As can be seen in Fig. 4, ( may - in principle 
- be as large as 0{2). However, for values of a^e'^^^'' ~ 0.5 x e~*^°°, which fit nicely 
into the picture developed in Subsection 3.3.1, we obtain (e'^i ~ 0.57 x e~*^°°, thereby 
suppressing pc to a sufficient extent. Interestingly, would then largely cancel 5 in 
(4.28). Moreover, also the first term in square brackets on the right-hand side of this 
equation may suppress pc even further. 



4.3 Refined Determination of 7 

Imposing the constraint that the values of Tc following from (4.25) and (4.26) agree with 
each other allows us also to refine the determination of 7 discussed in Subsection 3.4. 
In particular, it allows us to lift the degeneracy of the two solutions for 7 given in 
(3.64). In order to illustrate this feature, we show in Fig. 7 how the of ^ fit of the 
hadronic parameters to Acp{B(i — > tt+tt"), Acp{Bd — > n'^n') and H with and without 
a simultaneous fit to Tc varies as a function of 7. These fits are performed by fixing 7 and 
then fitting the hadronic parameters for each value of 7. We observe that the inclusion 
of the constraint from Tc does actually lift the degeneracy, leaving us with 

7 = (64.716:3)°. (4.29) 

The fact that the curve corresponding to the fit including Tc has its minimum in one of 
the minima of the other curve and practically touches the x-axis is non-trivial, but yet 
another sign of the consistency of our approach. 

Using Rij as given in (2.19) and the result for 7 in (4.29), the simple relations 

p — RbCOS^, fj — RbSin^, (4.30) 

which follow from the unitarity of the CKM matrix, allow us to calculate straightfor- 
wardly the other two angles a and P of the UT, where we obtain 

a = (93.61^°,=^)°, (3 = (21.711^)°. (4.31) 
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Figure 7: of a fit of tlie liadronic parameters to AQp{Bd — tt+tt^), JiQ^i^Bd — > tt+tt^) 
and H witli (daslied) and witliout (solid) a simultaneous fit to as a function of 7. By 
including the constraint from r^, the degeneracy is lifted. 

In Fig. 8, we follow [56], and compare these results with the allowed region for the apex 
of the UT following from the "standard analysis" [85].^ The solid window corresponds 
to the range for 7 in (4.29), whereas the dashed window indicates how the results change 
when the new Belle data (see Appendix B) are used. It should be noted that we show la 
regions, while the elliptic areas in the original UT plot show 95% CL. contours. Needless 
to note, the consistency of the overall picture is very remarkable. 

4.4 Analysis of the B — > ttK Observables 

Having all relevant hadronic parameters at hand, we may now analyse the observables 
provided by the B ttK modes within the SM, and may search for discrepancies that 
may shed light on NP. 

4.4.1 The Bd tt^K"^, iz^K System 

Let us first turn to the decays Bd ti^K^ and B^ —>■ tt^K. The characteristic feature 
of these modes is that EW penguins may only contribute in colour-suppressed form, as 
discussed in detail in Appendix D. Consequently, EW penguins are expected to have a 
marginal impact on the Bd tt^K^, B^ tt^K system. The relevant observables are 
the ratio R introduced in (1.5), which involves the CP-averaged branching ratios, and 
the corresponding direct CP asymmetries. If we assume again that pc = 0, (4.1) and 
(4.3) imply the following well-known expression [15]: 

= 1 - 2r cos5cos7 + r^. (4.32) 

^The small and large ellipses in Fig. 8 correspond to the SM and MFV, respectively, as obtained in 
an update [3] of [85]. 
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Figure 8: Comparison of our determination of 7 from the B tttt, hK data with the 
standard UT fit, where our la resuhs are indicated by the sohd and dashed windows, as 
described in the text. 



Using the SM value of 7 in (3.31) and the hadronic parameters given in (4.22), we obtain 



R\ 



SM 



n q4Q+0.033 



(4.33) 



which agrees well with the experimental result in (1.5). 

Additional information is provided by the — > n'^K^, B^ —>■ tt^K system through 
direct CP violation. As we have seen in (4.11), the CP asymmetry Aqp{B^ — > tt^K) 
vanishes for pc = 0, in accordance with the average of the current i?-factory data given 
in (4.12). Consequently, we are left with the direct CP asymmetry of the Bd tt^K^ 
modes, which we addressed already in (3.65) and (3.66). In terms of r and 6, this 
observable is given by 



•Acp{Bd 



TT 



2r sin 6 sin 7 



1 — 2r cos 6 cos 7 + 



(4.34) 



If we employ again - as in (4.33) - the SM expectation for 7 in (3.31) and the values of 
r and S in (4.22), this expression yields the prediction 



Ism 



n 140+0-139 



(4.35) 



which is in accordance with the current 5-factory average given in (3.65). Although 
we find that re**^ is strongly affected by non-factorizable effects, which is in particular 
reflected by the sizeable strong phase S in (4.22), the rather small value of r suppresses 
the direct B^ — > n^K^ CP asymmetry, in agreement with the data. 

To conclude our analysis of the B^ — >■ tt'^K^, B^ —>■ tt^K system, we emphasize 
that the corresponding i?-factory measurements give a current picture that does not 
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show any anomalous behaviour and is nicely consistent with the SM description. From 
this feature, we may in particular conclude that no anomalous enhancement of colour- 
suppressed EW penguins is indicated by the data. Further strategies to address these 
topologies are discussed in Appendix D. 



4.4.2 The Charged and Neutral B nK Systems 



Let us now analyse the data provided by the charged and neutral B — > tiK decays, where 
the quantities Rc and introduced in (1.3) and (1.4), respectively, play a key role. In 
analogy to R, they involve CP-averaged branching ratios. Using (4.1)-(4.4) with = 0, 
we obtain 



whereas 



with 



Rc — 1 — 2rc cos Sc cos 7 + r^, 

+qrc [2 {cos(5c + ^) cos — Tc coso; cos(7 — (/))} + qvc] , 



i?n = T [1 ~ 2r cos S cos 7 -I- r^l , 
b '- 



(4.36) 
(4.37) 



b — 1 — 2qrc cos(5c -|- to) cos -|- q 

+2p^ [cos^n COS 7 — qvc cos{9n — Sc — uj) cos(7 — 
The quantity b was introduced in [16] through 



(4.38) 



R^ 
R„ 



BR(^° ^ 7r°X°) + BR(^° ^ 7r°^°) 
BR{B+ TT+XO) + BR{B- TT-^o) 



rB+ 



1.19 ±0.16, (4.39) 



where the experimental value follows from [11]. This variable coincides with Rqq in [9]. 
In contrast to (4.32), (4.36) and (4.37) are significantly affected by the EW penguin 
parameter qe^'^e^^. Using the SM result 



qe'^'e" 



Ism 



= 0.69 X 



0.086 



(4.40) 



with iKife/Kfel given in (2.19) and the SM value of 7 in (3.31), the hadronic parameters 
in (4.22)-(4.24) yield 

(4.41) 



R. 



cISM 



1.14 



hO.08 
-0.07 



and 



R 



nISM 



1 1 1 +0.06 
-'-•-'--'--0.07' 



(4.42) 

exhibiting a pattern that is not in accordance with the current experimental picture 
given in (1.3) and (1.4), respectively. 

In this context, it is also interesting to consider the quantity 



(Rc - 1) + b(l - i?n) 

L = ^ ^-Z^ = 5.9 ± 2.8, 



2r2 



(4.43) 
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which was introduced in [16] and measures the violation of the Lipkin sum rule [86]; the 
numerical value in this expression corresponds to (4.26) and the averages of the S-factory 
data compiled in [11]. Using (4.36) and (4.37) with (4.38), we arrive at 

rlL = qrc [qvc - {tc coso; + pn cos(6'n - - '^)} cos(7 - 0)] 

+ [r cos 5 — Vc cos 4] cos 7 + Pn COS COS 7 + ^ (r^ — + p^) . (4.44) 

Taking, moreover, (4.17) into account, this expression can be simplified as follows: 
L = g [g — coso; cos(7 — 0)] 

Pn + rcos((5 - 6*11) 



— q cos(^n — — Lo) cos(7 



III analogy to (4.41) and (4.42), we may then calculate 

L|gM = 0.59l°:lt (4.46) 

which is in conflict with the experimental number in (4.43). 

Already back in 2000, when the observation of the — > 7r°ir° channel was an- 
nounced by the CLEO collaboration, two of us pointed out puzzhng features that were 
indicated by the corresponding results for and R^^ emphasizing the possibility of hav- 
ing NP in the EW penguin sector [14]. Now we have a much better experimental picture, 
where the BaBar, Belle and CLEO data show the same pattern for the observables R^ 
and Re- In a recent paper [16], we showed that enhanced EW penguins may in fact 
provide a solution to this puzzle. Here we perform a considerably more refined analysis 
and go beyond our previous study in the following respects: 

• We fix the relevant hadronic parameters re"', r^e!'^'' and PnC*^" through the 5 — > tttt 

analysis, as discussed above. In particular, we take also p„e*^" into account, which 
is found to be sizeable because of the very recent data for the B^, tt^tt'^ channel. 

• In addition to an enhancement of g, wc also allow for a CP-violating NP phase 
in the EW penguin sector, and for a CP-conserving strong phase a;, which could 
be induced by large non-factorizable SU (3)-breaking corrections. 

In order to make our analysis more transparent, let us first assume that c<j = 0°, 
as implied by the 5*^/(3) fiavour symmetry [13, 45, 48]. We are then left with the two 
EW penguin parameters q and 0. In Fig. 9, we follow [16] and consider the R^-R^ 
plane for the central values of the hadronic parameters in (4.22)-(4.24) and 7 = 65°, 
showing contours for different values of q and 6 [0°, 360°]. We observe that we may in 
fact move to the experimental region for an enhanced value of g ~ 1.8 and ~ —90°, 
where in particular the large CP-violating phase is in stark contrast to the SM picture 
characterized by (4.40). In order to put these observations on a more quantitative level, 
we use the value of 7 in (3.31) and the hadronic parameters in (4.22)-(4.24), which allow 
us to convert the experimental results for R^ and R^ in (1.3) and (1.4), respectively, into 
values of q and 0. Following these lines, we obtain 

5=1.751^:^9^ 0=-(85l}^)°, (4.47) 
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Figure 9: The situation in the Rn^Rc plane. We show contours for values of g = 0.69, 
q = 1.22 and q = 1.75, with G [0°, 360°]. The ranges from (1.3) and (1.4) (experiment) 
as well as (4.41) and (4.42) (SM) are indicated in grey. 

where our treatment of errors is discussed in Appendix C. We may now also calculate 
the quantity L in (4.45) within our NP scenario, yielding 



(4.48) 



in accordance with the experimental result given in (4.43). 

In addition to the CP-violating asymmetries that we considered already in our anal- 
ysis, there is yet another one that is strongly constrained by the i?-factory data [11]: 



BR(S+ ^ 'K^K+) - BR{B- TT^K' 



BR{B- 



TT' 



BR{B- 



TX 



0.00 ± 0.07. (4.49) 



Using (4.2) with pc = 0, we obtain 



Aqp{B — >• IT K ) = — [tc sin Sc sin 7 — qrc {sm{5c + uj) sin (f) + Vc sin 00 sin(7 — 0)}] , 

Rq 

(4.50) 

where the expression for R^ is given in (4.36). We are now in a position to calculate this 
CP asymmetry, where we obtain 



within the SM, and 



A'^^B^ - 



„0^±^| _QQi+0.14 

''Ism ~ '^•'^-'--0.10 



n'K^) = 0.04l°:i 



(4.51) 
(4.52) 



for our NP scenario, in accordance with the experimental value given in (4.49). On the 
other hand, the analysis performed in [16] favoured generically larger CP asymmetries. 
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Consequently, the determination of the hadronic parameters through the B — > tttt analy- 
sis and the introduction of the CP- violating NP phase 4> allow us now to achieve a much 
better agreement with the experimental picture. 

Interestingly, we may also employ the experimental information on Rc, Rn and the 
CP asymmetry Aqp{B^ n^K'^) to determine the three EW penguin parameters q, 
and uj simultaneously. This analysis gives 

q^l.7Atlil, 0=-(85ll^)°, cv^-iAttir, (4-53) 

where a central value of a; being nicely consistent with 0° - as expected in the strict SU (3) 
limit - is very remarkable. While this feature supports also our working assumption 
of using the SU{3) flavour symmetry, the reduction of the error on cu would be very 
desirable. 

Let us finally note, for completeness, that there is also a mathematical solution with 

- 0° and a; ~ 90°, 

q^2.32t\il = -(lOll^)°, c^ = (82li)°. (4.54) 

However, since such a value of to looks completely unrealistic, we will not consider this 
solution further. 



4.4.3 Elimination of the Second Solution for xe'^^ in (3.43) 

Let us now come back to the second solution for xe*^ represented by (3.43). Using (4.22) 
and the relations in (4.17)-(4.20), we obtain 



p^ = o.iit'oZ e^ = +{Q^tltr (4.55) 



re = 0.20^°:°?, = +(68li)°, (4.56) 

where the agreement between the values of Vc in (4.24) and (4.56) is obvious from (4.25). 
If we assume that u; vanishes, (1.3) and (1.4) give 

q=l-r3tlil cl>=-{102tlfr, (4.57) 

yielding 

A'^l>{B^ ^ n'K^) = 0.8ll°:°2^ (4.58) 

which is very much disfavoured by the experimental result given in (4.49). Even the 
inclusion of a non-zero strong phase cu does not provide any physically attractive solution 
for the EW penguin parameters, implying, for instance, values for q as high as 2.4 that 
are totally exluded by the rare-decay constraints discussed in Section 5. Consequently, 
we will not consider the second solution for xe^^ in (3.43) further. 
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4.5 Prediction of CP Violation in — ^ tt^I^s 

The decay vr^ii's is a transition into a final state with CP eigenvalue —1. Using the 

standard formalism for the calculation of the observables provided by the corresponding 
time-dependent CP asymmetry [2], as discussed for — > tt+tt" in Subsection 3.1, (4.4) 
yields 

2 r 

Acp{Bd n^Ks) = T grcsin(5c + uj) sin0 



-pn {sin sin 7 — qVc sin(^n — Sc — uj) sm{j — (p)} 



(4.59) 



A^^{Bd ^ 7r°i^s) = [sin (f>d - 2qr^ cos{S, + u) sin(0d + 0) + q^r^ sm{(j>d + 20) (4.60) 

+2pn {cos 6'n sin(0d + 7) - Q^c cos(6'n - (5c - a;) sin(0rf + 7 + 0)} + Pn sin(0d + 27) 

where the expression for b is given in (4.38). In the special case of pn = and = 0°, 
we arrive at the following, well-known results [59]: 

A^^iBa ^ n'^Ks) = (4.61) 

A^^{Ba ^ TT^iTs) = - sin 0, = ^^^^(5^ ^ J/z^i^s), (4.62) 

where the average of the current B-factory data for the mixing-induced B^ — >■ J/tpKs 
CP asymmetry is given in (1.8). Recently, the BaBar collaboration reported the first 
results for these observables [24] , which are given by 

A^^iB, ^ n'Ks) = +0.401°:^ ± 0.10 (4.63) 

A^^{B, ^ 7r°i^s) = -0.48;°;^^ ±0.11. (4.64) 

Moreover, there is also the following measurement of the direct CP asymmetry of the 
B^ n^K^ channel available [11]: 

Acp{B^ 7r°is:°) = -0.03 ± 0.36 ± 0.09, (4.65) 

which is supposed to agree with (4.63). Consequently, these experimental numbers are 
expected to change significantly in the future. It will also be very exciting to see the 
corrsponding first Belle results. 

Following the analysis performed in Subsection 4.4.2, we may predict the CP-violating 
Bd — > TT^i^s observables. In order to illustrate the dependence on the EW penguin 
parameters, we employ - in analogy to Fig. 9 - the central values of the hadronic pa- 
rameters in (4.22)-(4.24), u = 0°, 7 = 65°, and show in Fig. 10 the contours in the 
A^^{Bd 7i^Ks)-A'^^{Bd n^Ks) plane corresponding to various values of q with 
G [0°, 360°]. If wc take the uncertainties both of the parameters in (4.22)-(4.24) and of 
7 in (3.31) into account, and assume again that uj = 0°, the SM expression (4.40) yields 

A^l.{Bd ^ 7r°i^s)|sM - +0-12lS:}3 (4-66) 
A^^iBd ^ ttO^s) = -0.86in?, (4.67) 
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Figure 10: The situation in the A'S'piBd 7i^Ks)-A'^'p{Bd ir^Ks) plane. We show 
contours for values of g = 0.69 to g = 1.75 and (p G [0°,360°]. The grey area shows 
the lower left-hand corner of the BaBar la range, whilst the central value is outside the 
plotted area. The central values of (4.66), (4.67) and (4.68), (4.69) are indicated by the 
labels "SM" and "NP", respectively. 

whereas our NP scenario, which is characterized by the values of q and in (4.47), 
corresponds to the following prediction: 

A'^UBd ^ Tr^i^s) = +0.05l°:2^ (4.68) 

A^^iBd ^ n'Ks) = -0.99^ni- (4-69) 

The measurement of these CP asymmetries will allow a crucial test of our NP scenario. 
By the time solid experimental numbers are available, the uncertainties of the parameters 
entering (4.68) and (4.69) are expected to be significantly smaller, thereby leading to 
much more stringent predictions of CP violation in B^ n^Ks-^^ Constraints for these 
observables were also recently derived, within the framework of the SM, in [87]. 

4.6 Future Avenues Offered by Bs Decays 

The physics potential of i?s-meson decays, which can be exploited at hadronic 5-decay 
experiments, i.e. at run II of the Tevatron [71] and later on at the LHC [72], provides 
interesting strategies to explore CP violation and to obtain insights into hadronic physics. 
These results will nicely complement the B^ ^ t^t^i t^K methods proposed above. 

4.6.1 Bs K+K- 

As pointed out in [53], the decay Bg —>■ K^K~ can be related to B^ tt+tt^ through the 
?7-spin flavour symmetry of strong interactions, thereby providing attractive simultane- 
ous determinations of 7, d and 6, as well as insights into fZ-spin-breaking effects. It will 

""^"It is evident from Fig. 10 that the error in (4.69) is accidentally small. 
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be very exciting to see whether we will arrive at a picture that is consistent with the one 
developed in Section 3. Interestingly, we may use the hadronic parameters determined 
there to make predictions for the CP- violating Bg — > K'^K~ observables with the help 
of the [/-spin flavour symmetry, which implies the following expressions: 

A^l,{Bg - K^K-) = 2e^singsin7 

^ e2 + 2ed cos ^ cos 7 + ^2 ^ ^ 

j^mbcm K+K-) = ^^^i'^s + 27) + 2ed cos 9 sin(0, + 7) + ci^ sin 0, ^ 

+ 2ed cos 9 cos 7 + ' 

where the B^-B^ mixing phase (pg is given by 

0, = -2A2^ = 2/3, (4.72) 

in the SM, with the numerical value of Ps in (2.20). Note that the expression for 
A'^p{Bs K+K-) agrees with the one for A'^p{Bd ir^K^) in (3.66). Using now 
the range for 7 in (3.31) and the hadronic parameters in (3.41), we obtain the following 
SM predictions: 

- ^^^~)|sM = 0-14tg:o9 (4-73) 

A^^iBg K^K-)\^^ = -O.lSt'o'ol (4.74) 

where the latter observable may be affected by NP contributions to B^-B^ mixing. By 
the time the CP-violating Bg K^K~ observables can be measured, more precise SM 
predictions will be available. 

Moreover, we may also explore the branching ratio of this channel. To this end, we 
use the quantity H introduced in (3.59). If we use again (3.31) and (3.41), we obtain 

^IsM = 7.011:^7, (4.75) 

which is in nice accordance with the numerical value in (3.61).^^ In order to be able to 
predict BR(i?s K^K~), we have to know the [/-spin-breaking factor \C'/C\, which is 
given - within the factorization approximation - as follows [53]: 



(4.76) 



In a recent analysis [83], this parameter has been calculated through QCD sum rules, 
with the following result: 

^' = 1.76tg:l^. (4.77) 

fact 



c 



If we now complement (3.59) with (4.75) and (4.77), and use the experimental result 
BR(Bd ^ TT+TT") = (4.6 ± 0.4) X 10-6 [11], we obtain 



BR(S, ^ K+K-) = (35l^g) X 10-^ (4.78) 



^^Because of the discussion in Subsection 3.4, this feature is of course not surprising. 
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The very large uncertainty reflects the fact that our analysis is rather insensitive to H 
(this can also be seen in Fig. 1), and therefore the predicted value (4.75) has a large 
uncertainty that propagates to (4.78). 

Alternatively, wc may - in the spirit of (3.61) - assume that the penguin annihila- 
tion and exchange topologies discussed in Subsection 3.3.2 play a minor role. Taking, 
moreover, factorizable 5'?7(3)-breaking corrections into account, we obtain 



BR(I?,, ^ A"+A'-) 



Mb, ^M^/Mb,,Mk/Mb,) Tbo 



X 



BsK 



{Ml- 0^ 



where 



Fs,.(M2;0+) 
Fb,k{MI-Q+) fMl 



Ml 



Ml 



Ml 



M2 



H 2 



Ml 

Bd 



Ml 
Ml 



— i.40_o.i4 



corresponds to (4.77) [83]. If we then use the experimental result BR(5d 
(18.2 ± 0.8) X 10"^ [11], we arrive at 

BR(S, ^ K^K-) = (35 ± 7) X 10"^ 



(4.79) 
(4.80) 

(4.81) 



A measurement of this branching ratio, which should soon be available from run II 
of the Tevatron, will be very interesting, allowing in particular valuable insights into 
(4.77). It is interesting to note that - in contrast to (4.78) - the predictions of the CP- 
violating observables in (4.73) and (4.74) are not affected by factorizable t/-spin-breaking 
corrections [53], i.e. do not involve a ratio of form factors as in (4.76). 



4.6.2 TT^KT 



The decay 
symmetry of strong interactions 



TT^K is related to the 



TT mode through the t/-spin flavour 
, which allows us to write 



1 + -re'^e''^ 
e 



(4.82) 



complementing the S° — > tt K'^ amphtude in (4.3). If we then combine the direct CP 
asymmetry 



2er sin (5 sin 7 



+ 2er cos 6 cos 7 



(4.83) 



with its Bd — >■ TT^K^ counterpart in (4.34), we may determine r and 5 for given values 
of 7."^^ Fixing then Tc through (4.26), we may use (4.17) to eliminate pnC*^" in the 
B^ n^K^ amplitude, so that Rn and the two Bd — > n^Ks CP asymmetries depend on 

^^If we consider, in addition, the ratio of the CP-averaged Bs tt^K^ and — > tt^K 
branching ratios, we may determine 7 as well [88]. This extraction involves, however, the ratio 
Fb,k{MI- Q+)/Fb^^{MI- 0+) of S'i7(3)-breaking form factors. 
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the strong phase 6c and the EW penguin parameters q, 4> (and a;). If we complement 
these observables with and the direct — > t^^K^ CP asymmetry, we may extract 
these parameters and may perform internal consistency checks. 

The advantage of this avenue is that it is not affected by the penguin annihilation 
and exchange topologies discussed in Subsection 3.3.2. Consequently, it will be very 
interesting to sec whether we will eventually arrive at a consistent overall picture. If 
we neglect the penguin annihilation and exchange topologies and use the SU{3) flavour 
symmetry, we obtain the simple relation 

A^^iBs ^ TT^K^) = A^^{B, ^ TT+TT-), (4.84) 

which should already provide important insights into this issue. 

Let us finally note that EW penguins enter the Bs modes considered above only in 
colour-suppressed form. In the case of = 0°, these topologies would not affect our 
Bg strategies at all, as becomes obvious from the discussion in Appendix D. On the 
other hand, for (f) ^ 0°, anomalously enhanced colour-suppressed EW penguins would 
manifest themselves in the corresponding data. Consequently, the Bg studies complement 
the strategies to address the colour-suppressed EW penguins discussed in Appendix D. 

4.7 Summary 

Before turning to rare K and B decays as well as e'/ e in the next section, let us summarize 
the main results of our analysis of the B —>■ ttK system: 

• Employing the SU{2) isospin flavour symmetry of strong interactions, we have 
given a paramctrization of the B — > nK amplitudes, which allows us to deal with 
CP- violating NP effects in the EW penguin sector. Moreover, the relevant hadronic 
parameters are introduced in such a manner that we may determine them with the 
help of the S — > tttt analysis performed in Section 3. To this end, we have to 
neglect penguin annihilation and exchange topologies, and have to employ the 
SU{3) flavour symmetry of strong interactions. 

• We find a remarkable agreement between the corresponding determination of Tc 
and the value following from an alternative strategy, i.e. we arrive at a consistent 
picture, which would be spoiled if the working assumptions specified in the previous 
item were not satisfied. On the other hand, if we impose the constraint that the two 
values of Tc agree with each other, we may refine the extraction of 7 discussed in 
Subsection 3.4. In particular, this additional input allows us to resolve the twofold 
ambiguity, leaving us with 7 = (64.7l;g;g)°, which is in excellent agreement with 
the SM expectation given in (3.31). 

• Having all relevant parameters at hand, we may analyse the B — > ttK observables 
within the SM. As far as the B^ — ^ n^K^, B^ n^K system is concerned, where 
EW penguins may only enter in colour-suppressed form and are hence expected 
to play a minor role, we arrive at a picture that is in nice agreement with the 
S-factory data and does not indicate any anomalous behaviour. On the other 
hand, the SM analysis of the observables Rc and R^ of the charged and neutral 
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B — > ttK systems, respectively, which are significantly affected by EW penguins, 
is in conflict with the pattern of the S-factory data. 

• Interestingly, we can resolve this "puzzle" with the help of NP in the EW penguin 
sector. In particular, we arrive at the values of g = l.TSlgigg ^^^^1 (p = — (85l{4)°, 
i.e. at enhanced EW penguins with a large CP-violating NP phase. Moreover, 
we flnd the strong phase u> — —(4^32)°, in agreement with the implication of the 
SU{3) flavour symmetry that this phase vanishes. Let us emphasize that this 
picture corresponds to small direct CP violation in n'^K'^, as indicated by 
the S-factory data. Finally, we may predict the CP-violating observables of the 
Bd 7r°Xs channel, which will provide a crucial test of our NP scenario. 

• The impact of colour-suppressed EW penguins on our analysis has been addressed 
in Appendix D. The current i?-factory data do not indicate any enhancement of 
these topologies. 

• The decays Bg K^K~ and Bg n'^K^, which are very accessible at hadronic 
i?-decay experiments, complement the analysis of the Bu,d t^t^i t^K modes pro- 
posed above in a variety of ways. In particular the latter channel allows us to avoid 
assumptions about penguin annihilation and exchange topologies. 

5 Rare K and B Decays and e' /e 
5.1 Preliminaries 

As discussed in Section 2, the rare K and B decays are governed by the functions X, Y 
and Z, for which the parametrization in terms of |C| and 6c was given in (2.8). In order 
to obtain elegant expressions for rare decays, it is useful to go one step further and to 
rewrite (2.8) as follows: 

X=|X|e^^^, y=|y|e^^^, Z^\Z\e'^^, (5.1) 

with 

\X\ = V5.52 g2 _ 0.42 g cos <j) + 0.01, 
|F| = V5.52 g2 _ 3.00 g cos + 0.41, 

\Z\ = V5-52 g2 _ 4.42 gcos (p + 0.88, 

It should be emphasized that only two independent parameters, q and 0, appear in these 
expressions. 

Next we deflne the following weak phases: 

/3x = /5-/3.-^x, Py = P-Ps-Oy, Pz = P-Ps-Oz. (5.5) 



tan6'x ~ 



tan 6' 



Y 



tan 02 



2.35 q sin (j) 
2.35^cos0-O.O9 

2.35 (7 sin p 
2.35gcos0-O.64 

2.35 qsmc/) 
2.35gcos(/)-0.94' 



(5.2) 
(5.3) 
(5.4) 
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In [12], we have suppressed (3s in all formulae for rare decays but we have included its 
effect in the numerical analysis. 

Finally, the numerical constants in the formulae below correspond to [89] and 

sin^ = 0.231, a = A = 0.224, (5.6) 

128 

with the first two given in the MS scheme. They are the same as in [39]. The values of 
the remaining CKM parameters have been given in Section 2. 

With the result for (g, 0) in (4.47) at hand, we can calculate the functions X, Y and 
Z and the weak phases Px, Py and Pz- Setting iKfe/Kbl = 0.086, that is q — q, we find 

\X\ ^ \Y\ ^ \Z\ ^ A.Stli . (5-7) 

While the central value of |X| is still compatible with the data on K ttuu and B — > 
Xg^d^iy, the central value of |y| violates the upper bound |y| < 2.2 following from the 
BaBar and Belle data on B ^ XsH~^n~ [90], and the upper bound on BR(i^rL — > 7T^e'^e~) 
in (5.63) from KTeV [91]. In addition, \Z\ is too large to be consistent with the data on 
e'/e, even if the hadronic uncertainties in this ratio are large. 

On the other hand, the B — > ttK data seem to signal the possibility of enhanced 
values of \Y\ and \Z\ and of large weak phases 6i. Consequently, we may still 
encounter significant deviations from the SM predictions for rare decays, while being 
consistent with all experimental data. In order to illustrate this exciting feature, we 
consider only the subset of those values of {q, 0) in (4.47) that satisfy the constraint of 
|F| = 2.2. Using (5.3) and expressing q in terms of |y| and 0, and subsequently varying 
(f) in the full range given in (4.47), we obtain 

q = 0.92t'oZ = -(85111)°. (5-8) 

These values are compatible with all data on rare decays and also with the B nK 
data. In particular, we find the values for R^, R^, AcU^^ ^ 7r°ii'±), AcU^d ^ tt^Kq), 
AQp{Bd n^Ks), L and b given in the third column of Table 2. To this end we have 
varied (j) in the range given in (5.8), keeping |y| = 2.2. We compare the result of this 
exercise with the values for the B — > nK observables obtained without the \Y\ — 2.2 
constraint and with the data. 

Proceeding in the same manner, we also find 

|C| =2.24 ±0.04, = -(105 ± 12)° (5.9) 

|X| = 2.17 ±0.12, = -(86± 12)°, /3x = (HI ± 12)°, (5.10) 

|y| = 2.2 (input), = -(100± 12)°, /5y = (124 ± 12)°, (5.11) 

|Z| =2.27 ±0.06, = -(108± 12)°, /3z = (132 ± 12)°, (5.12) 

to be compared with C = 0.79, X = 1.53, Y = 0.98 and Z = 0.68 in the SM for 
rrit = 167 GeV. We will now turn to the implications of these results for rare K and B 
decays. 
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Ouantitv 


Without RD constraint 


With RD constraint 


Experiment 


Rr 


1 17 ± 12 (inuut) 




1.17 ± 0.12 


J. tjjl 


0.76 ± 0.10 (input) 


— 0.11 


0.76 ± 0.10 




'-'■'-'^-0.28 


o.o3;[J:i 


0.00 ±0.07 


A^S^iB.^TT^Ks) 


U.UO_Q 29 


0.08^°i^ 


0.40^^3^ 


A^^iBd^TT^Ks) 


_n qq+0.04 


-0 98+°-°^ 


-0.48l°:l« 


L 


fi 02+'^-'^'^ 

"•"-'^-4.67 


2 67+0-34 

^•^ ' -0.30 




b 


1 24+°-^^ 

-■-■^^-0.15 


1 ii;+0.16 
-'^■-'^^-0.13 


1.19 ±0.16 



Table 2: Theoretical predictions with and without the ll"] = 2.2 rare decays (RD) constraint 
and experimental data for the most important observables. 

5.2 K — > niyp 
5.2.1 Basic Formulae 

The rare decays tt+z/z/ and Ki^ tt'^uu proceed through Z^-penguin and box 

diagrams. As the required hadronic matrix elements can be extracted from the lead- 
ing semileptonic decays and other long-distance contributions turn out to be negligible 
[92], the relevant branching ratios can be computed to an exceptionally high degree of 
precision [93]-[95]. In [31]-[33], these decays have already been discussed in the NP 
scenario considered here. Below we update the formulae of these papers, adapt them to 
our notation, derive a few new ones, and include C(A^) terms that were neglected there. 
The branching ratios are then given as follows: 

BR{K+ TT+iyj?) = 4.78 X 10"^^ x [f^A'^Rf\X\'^ + 2fP^{X)A'^Rt\X\ cospx + PciX)^] 

(5.13) 

BR(Xl ^ Tr%P) = 2.09 X 10"^° x f^A^i?^^ |X | ^ sin^ (5.14) 



P,(X) = (1-^1 P,{X), (5.15) 



with 

PJX) = ^1 - 2 

where PdX) = 0.39 ± 0.06 results from the internal charm contribution [93, 94], Px is 
defined in (5.5) and f in (2.16). 

Once BR(i^+ -k^vv) and BR(i^^L t^^vv) have been measured, the parameters 
|X| and (5x can be determined, subject to ambiguities that can be resolved by considering 
other processes, such as the non-leptonic B decays discussed before and the rare decays 
discussed below. Combining (5.13) and (5.14), we find 

eWB, -B2- PcjX) 
= = ^cotpx: (5.16) 

£2V-D2 

where = ±1, and consequently 

2r 

sin2/?x = — ^. (5.17) 
l + ri 



s 
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Moreover, 



1^1- ^ £2sin/?x>0. (5.18) 



The "reduced" branching ratios Bi are given by 

_ BR(ir+ ^ T^^vv) _ BR(irL - ttV^) 

4.78x10-11 ' 2.09x10-10 " ^^"^^^ 

The formulae (5.17) and (5.18) were aheady presented in [32]. They are vahd for 
arbitrary 7^ 0° and generahze the ones given for the SM and MFV models in [96] 
and [97], respectively. For Qx — 0° and — — 1, one obtains from (5.17) the SM 
result that differs at first sight from the one given in [96]. However, in that paper and 
subsequent studies in the literature, a formula for sin 2/3 and not sin2(/9 — /9s) was given. 
Using the fact that = one can verify that (5.17), while being slightly more 

accurate, is numerically very close to the formula of [96]. In the scenario considered here, 
we have 99° < I3x < 125° and, consequently, Si — —1 and £2 = 1- 

It should be stressed that sin2/3x determined this way depends only on two mea- 
surable branching ratios and on Pc{X), which is completely calculable in perturbation 
theory. Consequently, this determination is free from any hadronic uncertainties and its 
accuracy can be estimated with a high degree of confidence. As in our scenario P and Ps 
are already known from the usual analysis of the UT, the measurement of in X — > tti/P 
decays will provide a theoretically clean determination of 9x- Similarly, a clean deter- 
mination of \X\ is possible by means of (5.18), so that (5.2) will allow us to determine 
q and (p. Assuming that the measurements of BR(fr+ n^vv) and BR(i^^L ti^vv) 
with 10% accuracy will be performed one day, the decays in question will most prob- 
ably provide the cleanest measurements of q and 0. This determination could then be 
compared with the one from other decays, in particular from AQl>{Bd — > n^K^) and 
A^^{Bd T^^Ks) that we proposed in Section 4. It could also be used to calculate 
other B ttK observables. 



5.2.2 Numerical Results 

Using the results for q and 4> in (5.8) and the parameters in (2.19), (2.20) and (5.6), we 
find 

BR{K+ TT+i/p) = (7.5±2.1) X 10-^\ BR(Xl ^ 7r°i/p) = (3.1±1.0) x 10"^°. (5.20) 

This should be compared with the SM prediction, for which we find 

BR(i^+ ^ 7r+i/i/)sM = (8.0 ± 1.1) X W-^\ BR{Ki^ 7r°z/i/)sM = (3.2 ± 0.6) x 10'^^ 

(5.21) 

in the ballpark of other recent estimates [92, 98]. On the experimental side, the results 
of the AGS E787 [99] and KTeV [100] collaborations are 

BR{K+ TT+uu) = (15.71^^2^) X 10"^^ and BR{Kj^ 7r%u) < 5.9 x 10"^ (5.22) 

respectively. 
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We observe that in our scenario BR(i^"'" ir^ui?) does not significantly differ from 
the SM estimate because the enhancement of the first term in (5.13) is to a large extent 
compensated by the suppression of the second term (cos/3x <S cos(/3 — Ps)) and its 
reversed sign. Consequently, BR(i^"'" tt'^vP) is very strongly dominated by the "top" 
contribution given by the function X. 

On the other hand, we observe a spectacular enhancement of BR(i^'L n^i'i') by 
one order of magnitude. Consequently, while BR(i^rL 7r°i/p) ^ (l/3)BR(i^'+ — > n'^uP) 
in the SM, it is substantially larger than BR{K'^ — > n'^uu) in our scenario. The huge 
enhancement of BR{Ki^ — > 7r°i/p) found here is mainly due to the large weak phase 
px ~ 111°, as 



BR(irL ^ vrVz/) 
BR(Kl ^ 7rVi/)sM 



X 



SM 



sin P 



x 



sin(/3 - Ps) 



(5.23) 



Inspecting (5.13) and (5.14), we observe that the very strong dominance of the "top" 
contribution in these expressions implies a simple approximate expression: 

fp^i;:;|.4.4x(.n&f.4.2±0.2. (6.24) 

We note that BR(ii'L — 7r°i/p) is then rather close to its model-independent upper 
bound [101] 

BK{Ki^ vrVz/) < 4.4BR(ii'+ ^ vr+t/p). (5.25) 

It is evident from (5.16) that this bound is reached when the reduced branching ratios 
Bi and B2 in (5.19) are equal to each other. 

A spectacular implication of these findings is a strong violation of the relation [96] 

(sin 2/5)^^^ = (sin2/3)^ii-s, (5.26) 

which is valid in the SM and any model with MFV. Indeed, we find 

(sin2/3),,p = sin2/3x = -(0.69^°^), (5.27) 

in striking disagreement with (sin2/3)^/^g = 0.736 ± 0.049 following from (1.8). 

In Fig. 11, we plot - in the spirit of [97] - BR(ii'+ tt+i/z/) as a function of BR(i^L 
n'^up) for fixed values of Px- As this plot is independent of \X\, it offers a direct 
measurement of the phase Px- The first line on the left represents the MFV models with 
Px = P — Ps, whereas the first line on the right corresponds to the model- independent 
Grossman-Nir bound [101] given in (5.25). The central value Px = 111° found here is 
very close to this bound. Note that the value of Px corresponding to this bound, where 
Bi — B2, depends on the actual value of these reduced branching ratios. As can be 
easily seen from (5.16), we have 

(cot/?x)Bound=-^^^^. (5.28) 
£2V B2 

For the central values of Pc{X) and B2 found here the bound corresponds to Px = 107.3°. 
As only cot Px and not Px is directly determined by the values of the branching ratios 
in question, the angle Px is determined only up to discrete ambiguities, seen already in 
Fig. 11. These ambiguities can be resolved by considering simultaneously other quantities 
discussed in our paper. The corresponding plot for different values of Px that are close 
to P can be found in [97]. 
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Figure 11: BR(K^ — n^ui?) as a function of BIi{Ki^ — > n^uu) for various values of 
The dotted horizontal lines indicate the lower part of the experimental range (5.22) and 
the grey area the SM prediction. We also show the bound in (5.25). 



5.3 B — ^ Xg and B 



s.d 



The inclusive decays B Xs^vv are also theoretically clean [93, 102]. Generalizing the 
known MFV formula to our NP scenario, we obtain 



BR(S ^ Xqvv) = 1.58 X 10" 



0.104 



0.54 

W). 



\V |2 

' \X(v)\\ 



(5.29) 



where q = s,d, and f{z) is a phase-space factor for B —>■ XcCv. The SM expectation 

BR(5 ^ XsUv)sM = (3.5 ± 0.5) x 10"^ (5.30) 

is to be compared with the 90% CL. ALEPH upper bound of 6.4 x 10~^. The exclusive 
channels are less clean but experimentally more easily accessible, with the 90% C.L. 
BaBar upper bound of 7.0 x 10^^. 

Next, the branching ratios for the rare decays Bq f^~^fJ'~ are given by 



BR{Bs fi^fi') = 2.42 X 10^^ x 



BR{Bd fi^fi') = 1.82 X 10"^ X 



1.46 ps 
1.54 ps_ 



238 MeV 



203 MeV 



\v;M'\yiv)\' 



\v:M'\Yiv)\' 



(5.31) 
(5.32) 



with tb^ ^ and Fb^ ^ being the lifetime and decay constants, respectively. 
In the SM, we find [103] 



BR{B, /x+/i")sM = (3.42 ± 0.53) 
BR{Bd /^'''/i~)sM 



AM, 



_18.0/ps 
1.00 ±0.14) X 10"^° 



X 10 



-9 



(5.33) 
(5.34) 
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This should be compared with the 90% C.L. bounds 

BR{Bs fJ.+fx-) < 9.5 (16) X 10"^ and BR{Bd < 1-6 x 10"^ (5.35) 

from the CDF (DO) and Belle collaborations, respectively [104, 105]. 

Inspecting formulae (5.2) and (5.3), we observe that for a fixed value of g, the effect 
of a non- vanishing phase 7^ 0° is to suppress \X\ and Thus, in spite of an 

enhancement of q, it is possible that for certain values of a suppression of the branching 
ratios given in (5.29), (5.31) and (5.32) relative to the SM expectations could be found. 
However, this does not happen in our case, and we find 



BR(B ^ X,vv] 



SM 



X 



SM 



2 



2.0 and BR(B,-,-V- 



Y 



SM 



2 



5.0, 



BR(B, ^ /x+At-)sM 

(5.36) 

with the same enhancements for B Xdi^v and B^ fJ-^fJ-^ ■ respectively. We also find 

BR{B X,vV) ^ 7xl0-^ BR{B, //V) ~ 17x10"^ BR{Bd //+//") ~ 5x10"^°, 

(5.37) 

which are still well below the experimental bounds. 
5.4 Kj, 

For the short-distance contribution to the dispersive part of Kj^ —>■ , we obtain in 
our NP scenario 

BR(Xl ^ A*+//")sD = 1.95 X 10"^ X [Pc{Y) +A'^Rt\Y{v)\ cos (3y]\ (5.38) 
where /3y is defined in (5.5), and 

^c(n=(l-y)n(n, (5.39) 

with Pc{Y) — 0.121±0.012 [94]. Unfortunately, because of long-distance contributions to 
the dispersive part of Ki^ fJ'~^fJ'~, the extraction of BR{Ki^ A^^Ai )sd from the data 
is subject to considerable uncertainties [106, 107]. While the chapter on this extraction 
is certainly not closed, let us quote the estimate of [106], which reads 

BR(Xl ^ /iV")sD < 2.5 X 10-^ (5.40) 

to be compared with BR(ii'L ^ /^+/^")sd = (0-8 ± 0-3) x lO"'^ in the SM. 

In the scenario with enhanced Z° penguins but no new weak phases considered by 

us in [16], BR(i^^L A*'''A*~)sd was substantially enhanced; in the present case, however, 
with 9y ^ 0°, the small value of cos/3y compensates the enhancement of so that we 
find 

BR{Ki^ /^V")sD = (0.9 ± 0.6) X 10"^ (5.41) 
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5.5 Forward— Backward Asymmetries in 6 — ^ sii^ii~ Modes 
5.5.1 Basic Formulae 

It has been pointed out in [34] that the exclusive decays Bd — >■ K*ix^ ix~ and their in- 
clusive counterparts B X^fi^fi^ offer excellent means to probe enhanced Z° penguins 
and, in particular, their complex weak phases. 

During the last two years significant progress for these transitions has been made both 
by experimentalists and theorists. On the experimental side, the Belle and BaBar col- 
laborations [90] reported the observation of this decay and of the B — > Xse'^e~ channel. 
The 90% C.L. ranges extracted from these papers [108] read as follows: 

3.5 X 10"^ < BR(B < 10.4 x 10^^ (5.42) 

2.8 X 10"^ < BR{B X^e+e") < 8.8 x 10"^ (5.43) 

On the theoretical side, important NNLO corrections have been calculated in [109]-[111]. 
The most recent reviews summarizing the theoretical status can be found in [112, 113]. 
The NNLO formulae are very complicated and it is not the purpose of our paper to 
present a detailed NNLO analysis here. For this reason, we give below NLO formulae 
that show transparently the size of various effects that we could expect in our scenario. 
Let us recall that these decays are dominated by the operators 

Qw = (s6)v'_A(/U/i)y, QioA = {sb)v-Aifi(J^)A, (5.44) 
which are generated through EW penguin diagrams. At low 

s = (Pl^^±pil, (5.45) 
ml 

also the magnetic operator Q77 plays a significant role. 

At the NLO level [114, 115], the invariant dilepton mass spectrum in the inclusive 
decay B — > Xg^i^ is given by 



d/dsV{h ^ sjj,^ 11 ) a 



2 



r(6 — > ceu) An"^ 



Vt 



ts 



ch 



f{z)K{z) 



U{s), (5.46) 



where 



U{s) = il + 2s) (^)r + ICiop) +4(^1 + ^^ |Cjf^|2 + 12C(f^ ReCfis), (5.47) 

and Cq^{s) is a function of s that depends on the Wilson coefficient Cg and includes also 
contributions from four-quark operators. Cj^j^^ is the Wilson coefficient of the magnetic 
operator Qjj. Exphcit formulae can be found in [114, 115]. The Wilson coefficients Cg 
and Cio are given as follows: 

Cg{v)^Po + ^^-AZ{v) + PEE{v), Cio(^) = -^^, (5.48) 
sm 0^ sm (fyf 
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with Pq = 2.60 ± 0.25 in the NDR scheme, Pe = O{10 ^) and v denoting collectively 
the parameters involved. Cq and Ciq are defined by 

C^viv) = ^Cg{v), Cwa{v) = ^Cio{v). (5.49) 

Of particular interest is the forward-backward asymmetry in 5 — > XsH~^n~ [116]. It 
becomes non-zero only at the NLO level. In our scenario it is given - in this approxima- 
tion - as follows: 



Afb(s) = -3Re 



(5.50) 



where U{s) is given in (5.47). The expression for the corresponding asymmetry in the 
exclusive decay K*fj,'^fi~ can be found in [34]. Both asymmetries vanish at a 

certain s = sq [117], which is determined in the case of the inclusive decay considered 
here through 

soReCf{so)+2Ci'^'^ = 0. (5.51) 

As in our scenario Cg^{sQ) and Cj^^^^ are assumed to be marginally affected by NP 
contributions, the value of sq is bound to be very close to the one predicted in the SM 
[109, 110]: So = 0.162 ± 0.008 after the inclusion of NNLO corrections. On the other 
hand, as we will see below, the magnitude of A-psis) and its sign could be strongly 
affected by the size of Cio and its phase Oy- 

Equally interesting is the forward-backward CP asymmetry introduced for — > 
K*ii^li~ in [34] , with the explicit expression given in equation (63) of that paper. This 
asymmetry can be substantially different from zero above the cc threshold if Ciq contains 
a large weak phase. The integrated asymmetry is found then to be [34] 

I^AZ = / ds AZ{s) = (0.03 ± 0.01) X (5.52) 
^0.52 KeGio 

where the range of the integration is chosen to decrease the effect of non-perturbative 
uncertainties. 



5.5.2 Numerical Results 

Using the NNLO formulae of [113] and comparing with the Belle and BaBar data for 
the branching ratios in (5.42) and (5.43) allows us to obtain a conservative upper bound 

\Y\ < 2.2, (5.53) 

which we already used in our previous numerical analysis [16]. Taking this upper bound 
into account and the values of |y| and Oy in (5.11), we find the following results that 
were anticipated in [34] in the case of a large CP- violating weak phase in Ciq. Our 
scenario is a concrete realization of such a situation. The most important findings are 
the following: 
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As 

sgn(ReC'io) = — sgn(ReF) = — sgn(cos^y) (5.54) 

governs the sign of Ap^is), and Oy = 0° in the SM, in our scenario with Oy > 90°, 
the sign of Afb(s) is opposite to the one in the SM: the asymmetry is positive for 
s < sq and negative for s > so- The value of sq is essentially unaffected. 



• However, 



ReC 



10 





Y 




\Y 


SM 



COS 6\ 



-0.5, 



(5.55) 



(ReCio)sM 

which through (5.50) implies the suppression of Afb(s) in our scenario with respect 
to the SM case. 

• On the other hand. 



AA^l = (0.03 ± 0.01) X tan^y 



-0.06 ±0.02 
-0.11 ±0.04 
-0.34 ±0.11 



[Uy 
{0Y 

(Oy 



115°) 
105°) 
95°) 



(5.56) 



is still another spectacular effect of the large weak phase implied by the B — > ttK 
data. The large sensitivity of A.Ap^ to 9y for 9y = C(100°) offers a useful method 
for the determination of this phase. 



5.6 Kl 



The rare decay Ki, — > n e'^e" is dominated by CP-violating contributions. It has re- 
cently been reconsidered within the SM [118] in view of the most recent NA48 data 
on 7r°e+e~ and Ki^ — >■ 7r°77 [119], which allow a much better evaluation of the 

CP-conserving and indirectly (mixing) CP-violating contributions. The CP-conserving 
part is found to be below 3 x 10^^^. Moreover, in the SM the indirectly (mixing) CP- 
violating contribution and the interference of both CP- violating contributions dominate 
the branching ratio in question, while the directly CP-violating contribution alone is 
significantly smaller and in the ballpark of 4 x 10~^^. In our scenario, this pattern is 
significantly changed, the latter part becoming the dominant contribution. Indeed, sim- 
ilar to BR{Ki^ — > tt'^i'u), the directly CP- violating contribution to BR(i^L 
enhanced by more than one order of magnitude. 

Generalizing formula (33) in [118] to our scenario, we find 



TT^e'^e ) is 



BR(Xl 
where 



7r°e+e-)cpv = IQ-^^ x 



3 X 10-4 



± C'dir 



K, 



3 X 10-4 



ImAf 



= fA^X^'Rt, f 



Vt 



ts 



cb 



0.98, 



(5.57) 
(5.58) 



sin(/3 — (3. 

C„,ix = (15.7 ± 0.3)|a,|2, \a,\ = l.OSt^oH, (5.59) 
Qnt = 1.02 y^v VC~, Cdir = 0.56 (y?^ ± f,y). (5.60) 



50 



Here 



\Y(v)\ 

mv = [^'o + PeE{v)] sin(/3 - + ^^-^ sm/3y - 4\Z{v)\ sin/3z 



sin^ 9^ 

Y(v) 



Via 



sin 9„ 



■ sin /3; 



Y, 



(5.61) 



(5.62) 



where Pq = 2.89 ± 0.06 [120] and Pe is 0(10"^). The effect of the NP contributions is 
mainly felt in y^A., as the corresponding contributions in y-^v cancel each other to a large 
extent. 

The present experimental bound from KTeV [91], 



BR(/s:l ^ 7r°e+e") < 2.8 x 10"^° (90% C.L.), 
should be compared with the SM prediction [118], 

BR(i^L ^ 7r°e+e-)sM = (3.2^^:2) x 10"^^ 



(5.63) 



(5.64) 



The enhancement of sin jSy over sin(/9 — f3s) and of [^'('y)! over Y makes the directly 
CP- violating contribution the dominant part of the branching ratio, for which we find 

BR(Xl ^ 7r°e+e-)cpv = (7.8 ± 1.6) x 10'^^ (5.65) 
which is lower than the upper bound in (5.63) by only a factor of 3. 

5.7 e'/e 

The formula for the CP- violating ratio e'/e of [121] generalizes for the NP scenario 
considered here as follows: 



rA'RtX" X F,,{v), 



(5.66) 



with 



FAv) = [Po + PEE{v)]sm{P-/3s) + Px\X{v)\smPx+PY\Y{v)\sm/3Y + Pz\Z{v)\sm/3z. 

(5.67) 

Here the Pi encode the information about the physics at scales /i < 0{mti Mw), and are 
given in terms of the hadronic parameters 



p - r(1/2) 



121 Me V 



ms{mc) +md{mc) 



(3/2) 



121 MeV 



ms{mc) + md{mc) 



as follows: 



(5.68) 



(5.69) 



The coefficients rf\ rf' and rf^ comprise information on the Wilson-coefficient func- 
tions of the A^" = 1 weak effective Hamiltonian at the next-to- leading order [35] ; their 
numerical values for different choices of at a — rUc in the NDR renormalization 

MS ^ 



scheme can be found in [121]. The hadronic parameters B^^"^^ and B^"^'^' represent the 



j(3/2) 
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matrix elements of the dominant QCD penguin operator Qq and the dominant EW pen- 
guin operator Qg, respectively. The numerical values of the Pj are sensitive functions of 
Rq and Rs, as well as of as- 

On the experimental side, the world average based on the latest results from NA48 
[122] and KTeV [123] and the previous results from NA31 and E731 reads 

e'/e = (16.6 ± 1.6) x 10"^ (5.70) 

While several analyses made in recent years within the SM found results that are com- 
patible with (5.70), it is fair to say that the large hadronic uncertainties in the coefficients 
Pi still allow for sizeable NP contributions. The relevant list of references can be found 
in [121]. 

In the SM, Ox = 9y = Oz = 0° and therefore (3x = (3y = Pz = P - Ps- We then find 
[121] that with 

Po = 19.5, Px = 0.6, Py = 0.5, Pz = -12.4, Pe = -1.6, (5.71) 

corresponding to = 1.2, Rs = 1.0 and as{Mz) = 0.119, an agreement with the 
experimental data can be obtained. 

In our scenario, the first term in (5.67) involving Pq and originating dominantly 
in the matrix elements of the QCD penguin operator Qq does not contain any NP 
contributions, while the important negative last term involving Pz and being related to 
the EW penguin operator Qg is strongly enhanced. With the values in (5.71), a negative 
s'/s is then obtained. Thus for our scenario to be consistent with the data, the hadronic 
matrix element of Qq or equivalently Rq must be significantly enhanced over Pg. For 
instance for B^^^^ — 2.0, B^^^^ — 0.62 and ms(mc) + ms(mc) — 106 MeV, corresponding 
to Re = 2.6 and Pg = 0.81, one finds for as{Mz) = 0.121: 

Po = 48.0, Px = 0.65, Py = 0.73, Pz = -10.4, Pe = -4.9. (5.72) 

Setting Px, Py, Pz at their central vahics, we obtain 

e'/e = 15.2 x 10-^ (5.73) 

which is consistent with the experimental data but the result for e'/e is very sensitive to 
the actual values of the coefficients Pj and the angles Pi. As reviewed in [121], Pg = 0.8 
used here is consistent with the most recent estimates that give Pg = 1.0 ± 0.2. Values 
for Rq as high as needed here have been reported in [124, 125]. In particular the values 
for Pe = 2.2 ± 0.4 and Pg = 1.1 ± 0.3 found in [124], that within the SM would give 
e'/e substantially higher than the experimental data, would be very welcome within the 
scenario considered here. 

5.8 Summary 

In this section, we have demonstrated that the sizeably enhanced EW penguins with 
their large CP-violating NP phase - as implied by the B — >■ ttK data - have important 
implications for rare K and B decays, as well as for e'/e. We find several predictions 
that differ significantly from the SM expectations and could easily be identified once the 
data improve. The most interesting results of this study can be summarized as follows: 



52 



• An enhancement of BR(i^^L ir^ui?) by one order of magnitude without any 
significant change in BR{K^ — > n'^uiy), implying BR(i^rL — > n^uP) to be close to 
its absolute upper bound derived in [101]. 

• A spectacular violation of (sin 2/5)7ri/p = (sin2/5)^/^g [96], which is valid in the SM 
and any model with minimal flavour violation. 

• A large branching ratio BR(i^'L Tr^e'^e") = (7.8 ±1.6) x 10~^^, which is governed 
by direct CP violation in this scenario, as opposed to the SM, where indirect CP 
violation dominates [118]. 

• A strong enhancement of the integrated forward-backward CP asymmetry for 
Ba^K*ii+li-. 

• Enhancements of BR(i? Xg^d^u) and BIi{Bs^d ~^ by factors of 2 and 5, 
respectively. 

• As far as e'/e is concerned, the enhanced EW penguins with their large CP- 
violating NP phase, as suggested by the B ttK analysis, require a significant 
enhancement of the relevant hadronic matrix elements of the QCD penguin opera- 
tors with respect to the one of the EW penguin operator Qg to be consistent with 
the e'/e data. 



6 Other Prominent Non-Leptonic S-Meson Decays 

As the last clement of our analysis, we explore the implications of our NP scenario of 
enhanced EW penguins with a new CP-violating weak phase for other prominent non- 
leptonic i?-meson decays, which play another key role in the physics programme of the 
B factories. 



6.1 B ^ (j)K 

In the SM, the B (t)K system is governed by QCD penguins [21] and receives sizeable 
contributions from EW penguin topologies [22, 23]. Consequently, these modes offer an 
interesting tool to search for signals of NP, and may of course also be affected within our 
specific scenario. 



6.1.1 Observables 

In order to address this exciting issue, we follow [29], and consider first 

1 — ■A(l)K 



with 



A. 



<t>K 



BKiB- 



6K^ 



1 + A^k'' 

+ BR{B- 



6K- 



BR(BO ^ + BR(BO ^ 4)K^) 



(6.1) 



(6.2) 
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and 

T^-^K ^ \ {AtA^^i - ^Ks) - A^l>{B^ ^ cf>K^)] (6.3) 

^ I [-^cpiB, - <PKs) + A'^l>{B^ ^ <PK^)] . (6.4) 

If wc use the experimental results given in [11] and add the errors in quadrature, we 
obtain 

_ f -0.05 ± 0.10 (BaBar) , . 

_/ -0.17 ±0.20 (BaBar) 

\ +0.08 ± 0.16 (Belle) ^ ^ 

+ _ r -0.21 ±0.20 (BaBar) 

\ +0.07 ±0.16 (Belle), ^ 

whereas the S-factory averages correspond to 

B^K = 0.00 ± 0.08, V-j^ = -0.01 ± 0.13, Pj^ = -0.04 ± 0.13. (6.8) 

We note that these observables are consistent with zero and hence do not indicate any 
deviation from the SM picture. In particular, as discussed in detail in [29] , B^k and "D^^ 
arc sensitive to the 1=1 isospin sector, thereby indicating that the corresponding ampli- 
tude f 1 is in fact suppressed, as is expected on the basis of plausible general arguments. 
On the other hand, T>'^j^ is sensitive to NP in the 1 = isospin sector, which involves an 
amplitude Vq with the sine of a CP-conserving strong phase Aq. The same contribution 
governs also the difference of AQp{Ba — > (pKs) and AQp{Ba J/i/jKs), where it enters, 
however, with the cosine of Aq. Within the SM, this difference is expected to vanish to 
a good approximation [26]-[30], whereas the current 5-factory data in (1.7) and (1.8) 
yield 

A'S^iB, ^ cpKs) - A^^iB, ^ J/i^Ks) = | ^ ^^f^^^^ (6.9) 

Consequently, we may well arrive at a discrepancy with the SM in the future, although 

the experimental picture is unclear at the moment. Should such a discrepancy really 
emerge and T>^j^ still be consistent with zero - in addition to Bcj)K and T>^j^, wc would 
have an indication for NP in the / = isospin sector, with a small CP-conserving strong 
phase Aq relative to the SM contribution. 



6.1.2 NP Analysis 

Since our scenario of enhanced EW penguins with a CP- violating weak phase cf) belongs to 
this category of NP, let us explore its implications in more detail. To this end, we neglect 
the amplitude vi of the 1=1 isospin sector, which corresponds to B^j^x — T^^k — 0, and 
write [29] 

^(^0 ^ 0i?°) = Ao [1 + ^;oe^^^°-*)] = A{B- ^R-) (6.10) 
A{Bl ^ 0^°) = Ao [1 ± ^;oe'(^°+'^)] = A{B+ ^ 0K+), (6.11) 
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where the CP-conserving strong amphtude Aq describes the QCD penguin contribu- 
tions, Vq measures the strength of the EW penguins with respect to the QCD penguins, 
which are governed by the (ss) (bs) quark-flavour structures having / = 0, and Aq is the 
CP-conserving strong phase mentioned above. Applying now once more the standard 
formahsm for the evaluation of the CP- violating observables, (6.10) and (6.11) yield 

2vq sin Aq sin </> 



Adir I D± 



1 -F2'yocosAocos0-F'u^J ' 
sin 0d 2vq cos Aq sin(0d 0) sin(0d 20) 



(6.12) 



(6.13) 



1 -I- 2vo cos Aq cos + 

For a simple ordcr-of-magnitude estimate of voe*^° within the SM, where = 0°, we 
assume that is dominated by internal top-quark exchanges, and use naive factorization 
with the leading-order Wilson coefficients given in [126] for = 225 MeV. Following 
these lines, we arrive at 



foe 



SM 
fact 



2(C9(mb) + Cio(mb)) 



4(C3(m,) + Q(m,))- 
2 X (-1.280 + 0.328) 



3C5(mfe) + Ce{mb) 
1 



X 



Acp{Bd 
If we use cos 



— = -0.20, (6.14) 

4 X (0.014- 0.030) + 3 X 0.009- 0.038J 128 ' v y 

which corresponds to Ao|fact ~ 180°; for a more refined treatment addressing also QCD 
penguins with internal up- and charm-quark exchanges, see [22]. In our scenario of NP, 
the enhancement of Cg + Cio is the same as that of our EW penguin parameter q. Taking 
into account also the constraints from the rare-decay analysis in Section 5, we conclude 
that this enhancement may be at most ~ 1.3 with respect to the SM, corresponding to 
values of Vq at the 0.25 level. We may then expand the CP-violating observables in vq, 
which gives 

V+j^ = -2vosmAosm(f) + 0{v^) (6.15) 
(PKs) - A'S^iBd ^ J/ijKs) = -2i;ocos Aosin0cos0d + 0{vl). (6.16) 

I = 0.68, Vq ~ 0.25, which suffers from large hadronic uncertainties, and 
take into account that our B — > ttK analysis points towards values of around —90°, 
corresponding to sin f=i — 1, we obtain 

Pj^ 0.5 X sinAo (6.17) 

A^^iB, ^ (l)Ks) - A^SpiB, ^ J/^Ks) ^ 0.3 x cos Aq. (6.18) 

Consequently, the experimental value of T>^j^ in (6.8) favours a small value of sin Aq. If 
we assume that the sign of cos Aq is negative, as in factorization, we conclude that the 
difference in (6.18) is negative as well. Since in our notation both mixing- induced CP 
asymmetries would equal — sin 2/3 ^ —0.7 in the SM, we generically arrive at 

(sin2/?)^^3 > (sin2/?)^K3, (6.19) 

where (sin 2/5)^^3 around +1 may well be possible, whereas the current result of the 
Belle collaboration in (1.7) points to a value with the opposite sign around —1. From 
the theoretical point of view, there are the following two possible loopholes: 
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Figure 12: The contours in tlie plane of A^^{Bd (pKs) - A^^{Bd J/^pKs) and 
"DJ^ for various values of vq with Aq G [0°,360°] and = —85°, corresponding to the 
central value of (5.8). 

i) In the case of the unconventional solution 0^ = 133°, the sign of cos0d would be 
negative [55, 56]. However, as we have noted in Subsection 3.4, this possibihty 
appears to be disfavoured. 

ii) We may encounter large non-factorizable effects, although it does not look likely 
that they may flip the sign of cosAq, which would require Aq around 0°. On the 
other hand, for Aq around ±90°, the impact of NP on AQp{Bd (pKs) would be 
very small, whereas 'D^j^ would be large, in conflict with the data. 

It is useful to consider the plane of Acp{Bd 4>K^) — A^p{Bd J/i/jKs) and V^j^, as 
done in Fig. 12, where we show also the current experimental results of the BaBar and 
Belle collaborations. Applying (6.12) and (6.13), we may calculate contours for different 
values of Vq, where each point is parametrized by the strong phase Aq. As far as the 
NP phase is concerned, we use = —85°, which is the central value in (5.8). In the 
future, plots of this kind will allow us to read off the hadronic parameters vo and Aq for 
our NP scenario easily from the improved B (pK data. 

6.2 B J/t/jK 

As emphasized in [27], EW penguins may have a non- negligible impact on the B J/ipK 
system. Let us, therefore, also address these effects in a more quantitative manner. The 
starting point is a systematic search for signals of NP arising at the decay-amplitude 
level of these modes [127], which is - from a formal point of view - completely analogous 
to the B (j)K system discussed in Subsection 6.1. Using the data for the CP-averaged 
B J/ipK branching ratios and the direct CP asymmetries compiled in [11] and [89], 
respectively, we arrive at the following picture of the B J/ipK counterparts of (6.1), 
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(6.3) and (6.4): 



B^K = -0.04 ± 0.04, = +0.023 ± 0.025, = +0.030 ± 0.025. (6.20) 

Consequently, these three observables do not indicate any deviation from the SM. 

Within our scenario of NP, we may write the S — > J/ ipK amphtudes in the same form 
as (6.10) and (6.11), with expressions for the direct and mixing-induced CP asymmetries 
as given in (6.12) and (6.13), respectively. Using again the naive factorization approach 
to estimate the relevant EW penguin parameter, we obtain 

. iAox |SM Cdimb) + Cio(to6) - (1 - l/A^c)Cio(mj,) 



Cgirrib) + Cio(mb) 



ibK 
0-2 



-1.280 + 0.328 
(125 



X — fsi -0.03. (6.21) 
128 ^ ' 



It should be emphasized that this expression suffers from large hadronic uncertainties. 
In particular, it is very challenging to estimate the effects of QCD penguins to the whole 
amplitude, which we may consider to be "effectively" included in the phcnomenological 
B J/'ipK colour-suppression factor af^. Nevertheless, it is instructive to consider 
again an enhancement of Cg + Cio of 1.3 and a NP phase (j) around —90°, which yields 

V-^j, ^ 0.08 X sin(Ao)v,ir (6.22) 

AcpiBd J/^pKs) ^ - sin 2(3 + 0.05 x cos(Ao)v,;f , (6.23) 

which should be compared with (6.17) and (6.18). Consequently, we may encounter a 
NP correction to the determination of sin 2/3 from the golden mode Bd J/ipKs at the 
0.05 level in our scenario, which corresponds to a shift of at most ±2° in /?. Such small 
effects are still beyond the present experimental and theoretical accuracy but could be 
reinvestigated in the LHC era. Moreover the comparison with similar effects in AMg^d 
could, in the future, shed some light on the size of the parameter {vQe'^^°)^K- It should 
also be noted that this parameter is expected to be different for higher J/ip resonances, 
so that the NP effects may largely cancel in the averages over such modes, which are 
usually formed by the BaBar and Belle collaborations and are also at the basis of (1.8). 

In order to search for NP effects in the "golden" mode B^ — * J/i/jKs, decays of the 
kind B(i Dn^, Dp^, ... are particularly interesting. If the neutral D mesons are observed 
through their decays into CP eigenstates D±, these pure "tree" decays, which do not 
receive any penguin contributions, allow extremely clean determinations of sin 2/3 [128]. 
Consequently, a future tiny discrepancy with (sin2/?)^i^g could also provide valuable 
insights into our specific scenario of NP. 



6.3 Summary 

The main results of our analysis of the B — > (j)K and B — > J/t/jK systems can be 
summarized as follows: 
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• Within our NP scenario, we expect generically (sin 2/^)0^:3 > (sin2/3)^Xs; where 
(sin 2/3)0x3 ~ +1 T^^y well be possible. This pattern is qualitatively different from 
the current S-factory data, which are, however, not yet conclusive. On the other 
hand, a future confirmation of this pattern would be another signal of enhanced 
CP-violating EW penguins at work. 

• The measurement of sin 2/3 through the golden mode — > J/ipKs may receive 
NP corrections as large as 0.05, which are, however, affected by large - essentially 
unknown - hadronic effects. The current 5-factory data disfavour an enhancement 
of the relevant NP parameter through such effects. 

7 Conclusions 

In this paper, we have developed a strategy that allows us to study simultaneously non- 
leptonic two-body B decays and rare K and B decays within the SM and its simple 
extension in which the dominant NP effects come from modified Z'^-penguin contribu- 
tions with a new CP-violating weak phase. This simple scenario is parametrized by only 
two variables (g, cf)) that take the values q — 0.69 and = 0° in the SM. 

The aim of our study was not only to learn about possible NP effects in the processes 
in question but also to determine the angle 7 and to gain some insight into the hadron 
dynamics, which plays an important role in non-leptonic B decays and e'/e, is less 
important in decays such as B ^ Xgii^ ii~ 1 and is essentially irrelevant in theoretically 
clean processes such as K ^ ttuP. 

If we had at our disposal precise experimental data for all processes considered here, 
the most efficient strategy would be to choose decays that are theoretically clean and free 
from NP contributions to determine the angle 7, use then this value in clean rare decays 
sensitive to NP in order to determine (^,0), and finally use 7 and (g, 0) in processes 
sensitive to hadron dynamics with the aim to obtain insight into the latter. In order 
to implement this procedure, clean tree-level strategies for the determination of 7 (see 
[2, 72, 128] and references therein), the rare decays K nuP, and non-leptonic decays 
such as -B — s> TTTT and B ttK could be used, respectively. 

Unfortunately, at present, the first two steps of this strategy cannot be made in 
view of the lack of data required for the clean strategies for 7 in question and in view 
of the insufficient experimental information on X — > nvP decays. On the other hand, 
data on S — > TTTT and B — > ttK are already available. Even though they are not yet very 
accurate, the progress expected in the coming years in measuring the relevant observables 
at the B factories allows us to expect that these non-leptonic decays will be known at 
an acceptable precision well before K — iruu decays - in particular Ki^ — tt^uu - and 
clean tree-level strategies for 7 will become useful. 

In view of this situation, it is necessary to concentrate first on the usual UT fits for 
7, S TTTT and B — > nK decays, determine 7, q, and hadronic parameters by using 
appropriate observables and SU (3) fiavour-symmetry and plausible dynamical assump- 
tions, and subsequently make predictions for rare decays. The three-step procedure for 
achieving this goal was presented in Section 1 and executed in the rest of the paper. 



58 



Our analysis has been summarized briefly in [12] and presented in detail in Sections 3- 
6 of the present work. A hst of results can be found at the end of each section. Here we 
recall the most important findings of our "anatomy" : 



• Using the available data for i? — > tttt decays and neglecting small EW penguin 
contributions to the relevant decay amplitudes, we were able to determine four 
hadronic parameters using only the isospin symiiietry of strong interactions and 
the information for the angle 7 from the UT fits. The result is given in (3.41) and 
(3.42). This determination is essentially free from theoretical uncertainties and the 
large errors in (3.41) and (3.42), which reflect the large experimental uncertainties 
in i?Jo^, A^I>{Bd TT+TT-) and A^^(Bd tt+tt"), should be significantly 
reduced in the coming years. Contour plots in Figs. 1 and 2 will allow a monitoring 
of these improvements. 

• Having the determined hadronic parameters of the S — > tttt system at hand, we 
may predict the CP-violating observables of the Bd — >■ 7r°7r° channel as given in 
(3.71) and (3.72), with the interesting perspective of having large direct and mixing- 
induced CP violation in this decay. Moreover, if at least one of these observables 
is measured, 7 can be determined in a theoretically clean way. 

• If we make the plausible assumption that penguin annihilation and exchange 
topologies play a negligible role and employ the SU (3) fiavour symmetry of strong 
interactions, we may complement the CP-violating observables of Bd — > tt+tt''" in 
a variety of ways with the data provided by the Bd — > tt^K"^ modes, which are 
only insignificantly affected by EW penguin contributions. Following these lines, 
we may extract 7, as given in (3.64). This determination can be refined with the 
help of other B — > ttK observables, with the result in (4.29), which is in very good 
agreement with the UT fits, as can be seen in Fig. 8. The remarkably consistent 
overall picture of our analysis also indicates that non-factorizable SU (3)-breaking 
effects are moderate and that our other dynamical assumptions are justified. 

• Using the hadronic S — > tttt parameters determined in the first step of our strat- 
egy, we may fix their B — > nK counterparts with the help of the SU (3) fiavour 
symmetry and plausible dynamical assumptions (see previous item), allowing us to 
analyse the B ttK system. Interestingly, we find that the SM cannot properly 
describe those B —>■ ttK observables that are sensitive to EW penguin contribu- 
tions. This is in particular the case for Rc and for which we obtain Rc ^ I.IA 
and i?n ~ 1-11 in the SM, whereas experiments give Rc ~ 1.17 and R^ ~ 0.76. On 
the other hand, the pattern of the 5-factory data for those B — > ttK observables 
that are only insignificantly affected by EW penguins does not show any anomalous 
behaviour, i.e. is in accordance with the SM picture. 

• We have demonstrated that all features of the present B irK data can be 
described within a NP scenario, where EW penguin topologies are moderatly en- 
hanced and carry a large CP- violating NP phase 4> around —90°, thereby requiring 
new sources for CP violation that he beyond the KM mechanism. Of particular 
interest are our predictions for AQp{Bd ^ Kg) and AQp{Bd^7r°Ks), and the 
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suggestion to use these observables in the future for the determination of the pa- 
rameters of the EW penguin sector. Moreover, also studies of Bg — > K'^K~ and 
Bs — > TT^K^ decays, which are very accessible at LHCb, complement our strategy 
nicely. We made also predictions for the corresponding S^-decay observables. 

• Restricting ourselves to a more specific NP scenario with enhanced penguins, 
we may explore the implications of our B — > ttK analysis for rare K and B decays, 
where in particular the CP-violating NP phase (f) has important consequences. In 
turn, the currently available rare-decay data have already some impact on the al- 
lowed ranges for the B — > tiK observables as summarized in Table 2. A detailed 
analysis of these effects has been presented in Section 5. Possibly the most in- 
teresting effects are found in the K — > ttuP system, where the enhancement of 
BK(Ki^ — > tt'^uP) by one order of magnitude and a strong violation of the MFV re- 
lation (sin 2/3)7rjyp = (sin 2/9)^^-3 are very spectacular. Similarly to the asymmetries 
A'^'piBd^Ti^Ks) and A^^{Ba^7i°Ks), the branching ratios for K+ ti+uu and 
Ki^ — > 7r°i/P will allow a useful determination of the EW penguin parameters and 
the comparison of these two very different determinations will be a very important 
test of the NP scenario considered here. 

• Finally, as discussed in Section 6, we expect a drastic modification of the Belle 
result (sin2/3)^Ks ^ (sin 2/3) ^^^-g, in that we find (sin 2/3)^x3 ^ (sin 2/3)^^-3. Con- 
sequently, within the NP scenario considered here, we find {sm2f3)^Ks to be of 
the same magnitude as the central value found by the Belle collaboration but of 
opposite sign! 

As we have seen in our analysis, studies of i? — tttt and B nK decays are not 
only interesting in the context of the exploration of CP violation and the search for 
NP, but also to obtain valuable insights into hadron physics. Consequently, improved 
measurements of the corresponding observables are also very important in order to see 
whether the theoretical approaches like QCDF [10], PQCD [129] and SCET [130], in 
addition to their interesting theoretical structures, are also phenomenologically useful. 
Independently of the outcome of these measurements, the phenomenological strategy 
presented here will be very useful in correlating the experimental results for B — > tttt 
and B — > nK with those for rare K and B decays, Ba — > (j)Ks and e'/e. 

Assuming that future, more accurate B — > tttt, ttK data will not modify significantly 
the currently observed pattern in these decays, the scenario of enhanced penguins 
with a large NP phase will remain an attractive possibility. While the enhancement of 
BR{Ki^ — > 7r°z/z/) by one order of magnitude would be very welcome to our experimental 
colleagues and {sm2l3)j^^p < would be a very spectacular signal of NP, even more 
moderate departures of this sort from the SM and the MFV expectations could be easily 
identified in the very clean K — > nuP decays as clear signals of NP. 
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A Compendium 



A.l The B ^ nn System 



1 + 2a; cos A + 



TDTTTT 

-'^00 



1 '2(1 cos 9 cos ^ + cP 

(P + 2dx cos( A — 6) cos 7 + 



Acp{Bd 



TT TT 



2d cos 9 cos 7 + 

2d sin sin 7 



1 — 2d cos ^ cos 7 + 
sin(0rf + 27) — 2d cos 9 sin(0d + 7) + 0?^ sin 0^ 



1 — 2d cos ^ cos 7 + 
2dx sm{9 — A) sin 7 



/I mix/ r? ^ ^0_i 



A.2 The B 



d^ + 2dx cos{9 — A) cos 7 + 
o\ _ + ^dx cos{9 - A) sin(0d + 7) + sin(0d + 27) 



+ 2dx cos{9 — A) cos j + x'^ 

TzK System 

R—l — 2r cos 6 cos 7 + 



1 — 2rc cos (5c cos 7 + 

+5rc [2 {cos(5c + to) cos — Tc coso; cos(7 — 0)} + 

i?n = 7 [1 — 2r cos 5 cos 7 + r^l 




b = — — 1 — 2qrc cos(5c + cos + r; 

itn 



2„2 
c 



2 



1 — 2r cos 5 cos 7 + 
2pc sin 9c sin 7 
1 + 2pc cos cos 7 + 



(A.l) 
(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 

(A.7) 

(A.8) 
(A.9) 



+2pn [cos COS 7 — QTc C0s{9^ — 5c — Oj) C0s(7 — (f))] + (A. 10) 

2/' sill () sill 



(A.ll) 
(A.12) 



Rr. 



Tc sin (5c sin 7 — grc {sin((5c + a;) sin + Tc sina; sin(7 — (/))}] 

(A.13) 

2 r 



A^piBd 77^ Ks) = - [grcsin(5c + a;) sm(f) 



-Pn {sin 9n sin 7 - qrc sm{9n - 5^ - co) sin(7 - 0)}j (A. 14) 
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^CP (^d ^ 7r°i^s) = [sin 0<i - 2qr^ cos{5^ + u) sin(0d + 0) + gV^ sin(0d + 20)(A.15) 
+2pn {cos 6'n sm(0rf + 7) - grc cos(6'n - 5c - ^) sin(0d + 7 + 0)} + Pn + 27) 

A.3 The Bs K+K', tt^K^ System 

Using the [/-spin flavour symmetry of strong interactions, we may express the corre- 
sponding CP- violating observables in terms of the hadronic — tt+tt" and ti^K^ 
parameters as follows: 

- K^K-) = 2e^singsm7 

' e2 + 2erfcos^cos7 + rf2 ^ ^ 

.mix/^ ^ = sin(0s + 27) + 2ed cos 6 sin(0, + 7) + sin 0, ^ 

^ e2 + 2ecicos^cos7 + rf2 ^ " > 



2er sin 5 sin 7 



(A.18) 



-|- 2er cos 5 cos 7 -|- 

B Comment on the New Belle Bd tt+tt" Results 

During the final stages of this work, the Belle collaboration announced the following 
update of the results for the CP- violating observables of the Bd 7r"'"7r~ channel [67] : 

Acp(Bd TT+TT-) - -0.58 ± 0.15 ± 0.07 (B.l) 
Acp{Bd TT+TT-) = -hl.OO ± 0.21 ± 0.07. (B.2) 

Using these new data, the averages in (3.29) and (3.30) change to 

A^^{Bd vr+TT") = -0.42 ± 0.13 (B.3) 
Acp{Bd TT+TT") = +0.70 ± 0.19. (B.4) 

Whereas the global picture of the analysis presented in this paper is not affected by 
these new numbers, some numerical results change. In particular, instead of (3.64), the 
averages in (B.3) and (B.4) would correspond to the following smaller values of 7: 

7 = (38.6l?i)° V (55.6^3^:;)°, (B.5) 

which would also be in accordance with UT fits, although a bit on the lower side. The 
contour corresponding to the new values was added to Fig. 8, and the central values in 
(4.31) change to a = 103.3° and P = 21.1°. 
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C Error Treatment 



Unless otherwise stated (e.g. in the determination of 7 in Subsection 3.4), we treat the 
errors in the following way: all predicted quantities depend on the hadronic parameters 
d, 9, X and A, and on some other input parameters. To take into account the fact that 
dependences on the hadronic parameters can cancel out, we do not vary the hadronic 
parameters inside their error bands, but rather vary the input parameters that are used 
to obtain them, i.e. A%p{Bd tt+tt-), ^cp (-^d Tr+Tr"), and R^^. Other input 
parameters that contribute to the errors are 7 and 0^, as well as A and additional 
hadronic parameters appearing in the context of Section 5 (some quantities involve also 
Rb) . The individual errors associated with the uncertainty of a specific input parameter 
are found by varying the corresponding parameter within its Icr band, while keeping 
the other parameters fixed at their central values. These errors are then added up in 
quadrature to obtain the total error for each observable. 



D Colour-Suppressed EW Penguins in B ^ ttK 
D.l General Structure 

In order to discuss the EW penguin contributions to the decays S+ — > tt'^K^ and 
n^K^, which are usually referred to as "colour-suppressed", we have to look at 
the EW penguin operators, exhibiting the following generic flavour structure: 



rw - \ Wc) - {ss) - (bb) + {2{uu) - (dd)}] (bs). 



(D.l) 



Using then the well-known isospin decomposition 

{uu) = - [{uu) + {dd)] +- [{uu) — {dd)], 
2 V ^ ' 2 s / 

1=0 1=1 



(D.2) 



we may decompose the EW penguin operators into isospin singlet and triplet pieces as 
follows: 

^ \2{cx:)-{ss)-{bb) + ^{{uu) + {dd)} 



{bs) 



(D.3) 



pen\ 



[{uu) - {dd)] {bs). 



:d.4) 



If we then apply the isospin flavour symmetry of strong interactions, as discussed in 
detail in [62], we obtain^^ 



A{B-' ^ 7r+X°)Ew = +2 
A{B', ^ 7t-K+)ew = -I 



pC{0) _ pC(l) 
EW EW 



C(0) 



EW 



R 



C(l) 



EW 



— g-^EW 
9 

— g-'^EW' 



:d.5) 



:D.6) 



^■^For simplicity, we suppress in the following discussion the primes introduced in Section 4 to distin- 
guish the B — > ttK amplitudes from their B — > tttt counterparts. 
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with 



pC(l) 



^eT = -^A^^ ^C'(X+7r-| [4(c-c) - 2(.-.) - 2(66) + {{uu) + (dd)}] (6.)|S°) (D.7) 



^..^^2 (D.8) 

summarizing the isospin singlet and triplet pieces, respectively. Here the combination of 
the generic Wilson coefficient C with the four-quark operators denotes symbolically the 
sum over the relevant EW penguin operators. We observe that 



C(l) 
EW ' 



(D.9) 



where the expression for P^w (-^-8) agrees with the one for the colour-suppressed 
EW penguin amplitude introduced in [46]. 

It is instructive to consider the tree-diagram-like matrix elements entering (D.7) and 
(D.8), which read as follows: 



dC{0) 
EW 



^0(1) 
EW 



v 2 2 



{uu){bs] 



d/T 



{uu){bs)\B',)T, 



and imply 



[-^Ew] T 



EW 



C(l) 



EW 



(D.IO) 



(D.ll) 



(D.12) 



The EW penguin operator Qg with the largest Wilson coefficient has the colour struc- 
ture {uaUa){hi3Sij) and its matrix elements (D.IO) and (D.ll) are consequently colour- 
suppressed. On the other hand, Qio ~ (waM/3)(6/3Sa) has a significantly smaller Wilson 
coefficient. The coefficients of the (V — A) ® (V -\- A) EW penguin operators Q-j and 
Qs are even further suppressed, so that these operators can be neglected. Consequently, 
the hadronic matrix elements in (D.IO) and (D.ll) are colour-suppressed, and (D.12) 
represents the picture of the colour-suppressed EW penguins that is usually adopted in 
the hterature. 

In our analysis of the P — > tttt system in Section 3, we have seen that penguin 
topologies with internal up- and charm-quark exchanges, which correspond to matrix 
elements with penguin-like contractions of {uu)(hs) and {cc){hs) operators with u and 
c quarks running in the loops [62], play an important role. Consequently, the standard 
picture of colour-suppressed EW penguins in P — > 'kK could also be affected through 
these penguin topologies. Moreover, the parameters 



ace 



pC- 
-^EW _ 


pC- 


Pew 


pC- 1 pA 
EW ^ EW 


pC+ 


pC+ 


-^EW _ 


-^EW 


Pew 


pC+ , pA 
-^EW "T -^EW 



(D.13) 
(D.14) 
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where and P^^^ are the colour-allowed EW penguin amplitudes contributing to 

V2A{B+ 7r°X+) and y/2A{B'^ n^K^), respectively, may not be as small as naively 
expected. In (D.13) and (D.14), we have used that the isospin flavour symmetry implies 

(D.15) 



EW 



pC- I pA 
-^EW ^ -^1 



EW 



EW 



EW' 



where Pew enters the EW penguin parameter qe^'^e''^ . As in Subsection 2.3, we define 
the weak phases of and P^-^ as the one associated with the combination Cg + Cio 
of Wilson coefficients, which can be done approximately. 

In the remainder of this appendix, we discuss in detail the possible impact of the 
parameters in (D.13) and (D.14) on the B nK analysis performed in Section 4, and 
propose strategies to search for indications of their effects in the corresponding data. 



D.2 Generalization of the Decay Amplitudes 

The generalization of the decay amplitudes in (4.1)-(4.4), taking the colour-suppressed 
EW penguin topologies into account, is given by 



A(B+ ^ 7r+X°) = -P' 
V2A{B+ 7r°X+) = P' 
A(P^ ^ TT-K^) = P' 



2 
3' 



with 



= -P' 



2 
3' 



1 + p^e'^-e'^ - ( 1 - ^ace'^c ) qe^'^e'^r^e'^^ 



"EW'^ 



qe Tee ace 



(D.16) 
(D.17) 
(D.18) 
(D.19) 

(D.20) 

(D.21) 

Because of the isospin relation in (D.15), the terms proportional to dc and ac have to 
cancel in the sums of the A{B+ n+K^), V2A{B+ n°K+) and A{B^ tt-K+), 
V^A{B^ — > 7r°ir°) amplitudes, respectively. 



D.3 The Case of = 0° 

Let us first consider = 0°, which applies also to the case of the SM. Although the 
analysis of the current P-factory data performed in Section 4 favours a large phase of 
(j) ~ —90°, it is interesting and instructive to have a detailed look at the situation arising 
for (f) = 0°. In this case, we may straightforwardly absorb the terms in (D.18) and (D.19) 
proportional to ace*^*^ in the amplitude P', yielding 

A{B^a n-K+) - P' [1 - re'^e'^] (D.22) 
v^A(P° ^ 7r°X°) = -P' [1 + pr^e'^-e'^ - qe^r^e'^'^] . (D.23) 
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As far as (D.16) and (D.17) are concerned, we obtain amplitudes of the structure 



A{B 

with 



+ 



TT 



+ = -P' 



1 + Pee' 



6 ^EW ^ 



(D.24) 



n^K+) = P' [l + pce^^^e*^ - |e*^ - (^1 - a[!^e^^c'^ qe^^ r^e'^'^j , (D.25) 



.0(1) ivg' 
"ew e 



pC(l) 
EW 

-Pew 

qe'''r^e'^'^a'i^e'^c\ 



(D.26) 
(D.27) 



where Pew"* given in (D.8). It should be noted that P', re^^, PnC*'^", TcC**'^ and pcc'^'^ 
now contain also contributions from colour-suppressed EW penguins. The expressions 
in (D.22) and (D.23) are the counterparts of (3.74) and (3.75), respectively. 

We observe that the colour-suppressed EW penguins do not explicitly affect the 
analysis of the P — > tttt modes and the neutral B — > tiK decays. Because of this feature, 
P° — > T^~K'^ is actually the "natural" partner of B^ — > 7r"'"7r~ to deal with the famous 
penguin problem and not the B^ — > n^K^ channel, as is sometimes done in the literature. 
If we use the B ^ tttt data as described above to fix the hadronic P nK parameters, 
we may extract q and uj from and one of the CP-violating P^ t^^Kq, observables; 
the remaining CP asymmetry provides an important consistency check of the = 0° 
scenario, in particular to see whether this weak phase actually vanishes. Interestingly, 
enhanced colour-suppressed EW penguins would yield a difference between the values of 
Tc following from (4.25) and (4.26), which is not indicated by the data. 

If we fix the hadronic P — ^ tiK parameters through the P tttt system, the analysis 
of the charged B — > tiK modes may be affected both by the colour-suppressed EW 
penguins and by the parameter. Following the avenue described above, we may predict 
the corresponding observables, and may check whether we obtain agreement with the 
experimental picture. Alternatively, we may use the value of Tc following from (4.26), 
and may shift the terms in (D.24) and (D.25) proportional to ag^e^^c into P' . We may 
then also analyse the observables of the charged P nK system in a manner that is 
not affected by the colour-suppressed EW penguins, whereas Pc may still enter. 



D.4 The Case of 7^ 0° 

In the case of 7^ 0°, we could of course still absorb colour-suppressed EW penguin 
contributions in the amplitude P'. However, as these terms now contain a CP-violating 
weak phase, we have no longer the simple amphtude structure of P' — |P'|e"^^", where 
5p> is a CP-conserving strong phase. 

If wc look at (D.16) and (D.18), we see that the observable R and the direct CP 
asymmetries of the P^ n^K and B^ n^K^ modes would be significantly affected 
by the terms proportional to ace*^^; and ace*^^; should these coefficients not be small 
quantities, i.e. should the colour suppression of the EW penguins not be effective. How- 
ever, as discussed in Subsection 4.4.1, these observables do not indicate any anomalous 
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behaviour. Moreover, as we have already noted above, also the agreement between the 
values of Tc in (4.24) and (4.26) does not favour an enhancement of the colour-suppressed 
EW penguins. In this context, it is also important to note that the B — > nK data do not 
require a dramatic enhancement of the parameter q. Moreover, the analysis of the rare 
decays in Section 5 favours q to be even smaller than the central value given in (4.47). 

In order to deal further with the colour-suppressed EW penguins in the 7^ 0° case, 
it is useful to introduce 



i{u>) 



SO that 



1 _ tace'^c 
3 ^ 



1 

— ( 
2 



'C 



iA 



(1) 



qe" 



1 ane 

3 ^ 



with 



''C 



{q)e 



P 



,C(0) 



EW 
^EW 



(D.28) 

(D.29) 
(D.30) 

(D.31) 



in analogy to (D.26). Considering only tree-diagram-like matrix elements, as in (D.IO) 
and (D.ll), yields 



SL 
.(1) 



1 

3' 



(D.32) 



If we assume that also a^^ is suppressed with respect to a^' , we may identify the terms 



(1) 



proportional to qe^ in (D.17) and (D.19) simply with the "effective" EW penguin pa- 

erestingly, the strong phase A[!^ n 
although cu could still be tiny. If we complement then 



rameter {q)e^^^K Interestingly, the strong phase A[!^ may induce a sizeable value of (u), 



BR(Z?- 



BR(B- 



Rc 
T 



0.99 ±0.15 



(D.33) 



with two of the three CP- violating observables provided by the i? — > n K^, — > tt i^s 
modes, we may determine (g), {ui) and 0; the remaining third CP asymmetry can be 
predicted and allows a crucial consistency check. 

To conclude, let us emphasize once more that our analysis of the current B-factory 
data for the B nK modes performed in Section 4 points towards a consistent overall 
picture. In particular, R and the direct CP asymmetries of the B^ n^K^, B^ — > n'^K 
modes do not show any anomalous behaviour. Consequently, we have no experimental 
indications for an enhancement of the colour-suppressed EW penguins, which have actu- 
ally a very complicated internal structure. In the future, the strategies discussed in this 
appendix will allow us to explore these contributions in a more stringent manner with 
the help of improved experimental data. 
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